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We reconsider in some detail a construction allowing (Borel) convergence of an alternative per- 
turbative expansion, for specific physical quantities of asymptotically free models. The usual per- 
turbative expansions (with an explicit mass dependence) are transmuted into expansions in 1/F, 
where F ~ l/g(m) for m S> A while F ~ (m/A) a for m < A, A being the basic scale and a given 
by renormalization group coefficients. (Borel) convergence holds in a range of F which corresponds 
to reach unambiguously the strong coupling infrared regime near m — -> 0, which can define cer- 
tain "non-perturbative" quantities, such as the mass gap, from a resummation of this alternative 
expansion. Convergence properties can be further improved, when combined with S expansion (vari- 
ationally improved perturbation) methods. We illustrate these results by re-evaluating, from purely 
' perturbative informations, the O(N) Gross-Neveu model mass gap, known for arbitrary N from 

| exact S matrix results. Comparing different levels of approximations that can be defined within our 

framework, we find reasonable agreement with the exact result. 
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In many quantum field models, non-perturbative results may be obtained from the 1/N expansion El S|, which at 
, leading orders resums only a certain class of graphs of the original perturbation series in the coupling. In parallel, 
for a given perturbative series there are more direct efficient summation techniques, like the Borel || Q or Pade|i) 
fS| | methods, as well as generalizations combining the latter two|| [jj. Typically, the Borel method is useful even for 
. non Borel-summable perturbative expansions, as for instance in QCD, since it gives interesting informations on the 
incompleteness of the pure perturbation theory, and the necessary additional non-perturbative (power corrections) 
contributions to a given physical quantity. Also, a rather different modification of the usual perturbation theory, known 
O ' as delta-expansion (DE) or "variationally improved perturbation" (VIP)|| 0, is based on a reorganization of the 
interaction Lagrangian such that it depends on arbitrary adjustable parameters, to be fixed by some optimization 
1 prescription. In D = I field theories, the quantum mechanical anharmonic oscillator typically, DE-VIP is in fact 
equivalent to the "order-dependent mapping" (ODM) resummation method^, and optimization is equivalent to 
a rescaling of the adjustable oscillator mass with perturbative order, which can essentially suppress the factorial large 
order behaviour of ordinary perturbative coefficients. This appropriate rescaling of a trial mass parameter was proven 
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to give a rigorously convergent series) If 12| e.g. for the oscillator energy levels |13[ and related quantities. 

In the present paper we reconsider yet another alternative expansion, proposed some time ago JT^|— 1 16 , which 



is close to (and partly inspired by) the DE-VIP idea, but particularly suited to apply gcncrically to arbitrary 
higher dimensional (renormalizable) models. Rather similarly with the 1/N expansion, it starts with a specific 
approximation but can include, at least in principle, the full information from the Lagrangian in a systematical way. 
The basic construction exploits a physically motivated renormalization group (RG) "self-consistent mass" solution, 
which resums RG dependence to all orders (at least in specific renormalization schemes). Moreover it provides 
a non-perturbative information, encoded in the infrared properties of an implicit function F(m), which may be 
viewed as a generalization of the logarithm, and occurring as the exact solution of the above mentioned RG equation 
with the self-consistent mass boundary condition. At first RG order, F(fh) = ln(m/A) — ^41nF(m) is essentially 
the Lambert function]!^], where A depends in a simple way on the mass and coupling RG coefficients (m is the 
renormalization scale-invariant Lagrangian mass and A the basic RG scale). Unlike the logarithm, however, F(m) 
has a power expansion behaviour in (m/A) 1 /' 4 in the infrared, for |m| < to c < A, while it matches the ordinary 
perturbative effective coupling, F ~ ln(m/A) ~ l/<?(m), for the short distance m ^> A perturbative regime. As we 
shall argue, F thus defines a rigorous (analytic) bridge between the usual perturbative short distance regime, and 
the strongly coupled, non-perturbative, massless (chiral) limit corresponding to to < A, or to <C A. This transmutes 
the ordinary expansion (in the coupling g) of physical (on-shell) Green functions, depending explicitly on a single 
mass to, into a 1/F expansion, or equivalently a (mass) power expansion in (to/ A) 1 /' 4 for sufficiently small to. The 
main idea is to use those properties of F, in asymptotically free theories (AFT), to infer a non trivial mass gap 
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M(to — > 0) 7^ Qp9 in the deep infrared (strongly coupled) regime, equivalently here the massless limit m — > 0. More 
generally other physical quantities can be derived similarly, for instance the quark condensate (qq)(m — > 0), one of 
the order parameter of chiral symmetry breaking (XSB ) in QCD|l5|, from their known perturbative expression 
for m>A. At this stage, it is important to remark that our construction is not by itself a proof of dynamical (chiral) 
symmetry breaking, and applies indeed independently of whether chiral symmetry is (dynamically) broken or not: it 
introduces rather an explicit chiral symmetry breaking mass m in the Lagrangian, in such a way that the properties 
of F(fh) encode a non trivial M/A(m — > 0) ratio, smoothly extrapolated from the massive case M(rh 0). This is 
to be simply viewed as a generalization, for to 7^ 0, of dimensional transmutation ]fl8| . 

Unfortunately, such an extrapolation to to — > is well-defined as far as the pure RG dependence of the relevant 
physical quantities is concerned, while it turns out to be badly afflicted when considering for the latter their 
complete perturbative (non-RG-dependent) expansion coefficients with their expected behaviour at large orders. As 
is well-known, in most models the leading large order behaviour of purely perturbative coefficients exhibit same signs 
(thus non Borel summable) factorial divergences (the infrared renormalon singularities [Q, [lj]]). The standard and 
seemingly unavoidable interpretation is that it implies large perturbative ambiguities, e.g. of 0(A) for the (pole) mass 
gap, reflecting incompleteness of purely perturbative expansions and the necessity of adding non-perturbative power 
corrections. However, our alternative expansion can be smoothly extrapolated down to small, and even negative 
(or more generally complex) values of the expansion parameter 1/F, thanks again to the properties of F, in such a 
way that the corresponding perturbative series can be Borel summable p0|. More precidely, in the simplest situation 
(corresponding to the RG parameter particular value A = 1), there is one branch of F such that F < 0, which 
simply produces the required sign-alternation in the perturbative coefficients ~ F~ n . In this particularly simple case 
A = 1, the range where F < happens to correspond also to to < 0. This is not a problem in principle, since in 
relativistically invariant theories the absolute sign of the Lagrangian mass term is irrelevant to physical quantities, 
moreover in our context the physically relevant results are in the massless (chiral symmetric) limit anyway. Thus 
taking Re[F] < (Re[m] < 0) simply corresponds physically to reach a strongly coupled regime near the massless 
Lagrangian limit, but without the usual ambiguities from renormalons. More generally, i.e. for an arbitrary AFT in 
an arbitrary scheme, as we shall examine there exist (complex) branches of F near the relevant massless limit, which 
is sufficient to ensure Borel summability, the Borel singularities being moved away from the real axis. Those different 
branches of F may correspond either to to > or to m < 0. We shall argue that the actual physical result is indeed 
independent of the branch on which the massless limit is reached, though the (unphysical) perturbative expansions 
have obviously different Borel summation properties depending on the branch of F considered. 

Independently of these Borel convergence properties of the 1/F-series, we can also consider, in a second stage, 
an appropriate version of the (order-dependently rescaled) "variationally improved" perturbation (DE-VIP), to 
be performed on the series in 1/ F(rh/ A), essentially replacing the true physical mass by an arbitrary adjustable, 
trial mass parameter. This produces a renormalization scheme (RS) dependent factorial damping of the original 

12]. Now here, the damping appears 



jut can further improve [EOl the Borel 



perturbative coefficients at large orders, similarly to the oscillator case] 11 
insufficient to make the DE-VIP series readily convergent for arbitrary to, 
convergence properties, which are also obtained typically for Re[F] < 0. 

The basics of our construction was defined before p4[, and some of its phenomenological applications explored, to 
some extent, in QCD ]l5|, ^6|. However, those previous numerical results were based on rather ad hoc approxima- 
tions, either by constructing approximants only based on the lowest orders of the perturbative expansion, and by 
optimization with respect to the renormalization scale and/or scheme |p"6f . In particular, it ignored completely the 
above mentioned ambiguities due to the large order, factorial behaviour of perturbative expansions. We are thus 
mainly concerned in the present paper to provide a more concrete illustration of the formal Borel convergence results 
obtained in ^0|, by considering the mass gap of the 0(N) Gross- Neveu (GN) model[^l], ^2|, known for arbitrary N 
values from exact S matrix results]^ and Thermodynamic Bethe Ansatz|24| similarly applicable p5| in many other 
integrable 2-D models. Those exact results serve as a test of our method, which is based in contrast only on the 
perturbative information, thus applying a priori to any other (asymptotically free) renormalizable models, e.g. in 
four dimensions, for which there are obviously no exact S-matrix results available. Typical examples are 4-D gauged 
AFT with nf massless fermions, like QCD, where the expected [p6| SU(rif)L x SU(nf)n — > SU{nf)y breaking is 
characterized via non-perturbative order parameters, generalizing the role of the mass gap in simpler models of 
dynamical symmetry breaking p7|. 

The paper is organized as follows. In section 2, we recall the main steps of our construction. We define the new 
perturbative expansion, keeping as much as possible the discussion general for any AFT but with specific illustrations 
in the O(N) GN model. We also give additional formulas and some properties which had not been discussed in 
refs. (l4|-[jl6| or 20 . In section 3, we briefly recall the usual problems of infrared renormalon singularities, the 



resulting perturbative ambiguities, and how these ambiguities are usually removed by non-perturbative contributions, 
whenever the latter can be explicitly evaluated, e.g. at the next-to-leading \/N order in 2-D integrable models like 
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the GN. In section 4, we reexamine the Borel convergence properties [|20| obtained for F < in the simplest situation, 
providing also a more detailed analysis of some technical issues. Section 5 and 6 analyse how the variationally 
improved perturbation (DE-VIP) can further improve these Borel convergence properties, and in section 6 are also 
introduced other (non-linear) convenient generalizations of the simpler DE-VIP construction, which can lead to a 
directly convergent alternative series. Finally in section 7 we give numerical applications for the O(N) GN model, 
where we compare in some details different approximations and/or rcsummation methods which can be constructed 
order by order within our approach. Section 8 contains some conclusions, and a number of technical issues used in 
various parts of the paper are discussed in five appendices. 



II. RG SELF-CONSISTENT MASS AND ALTERNATIVE EXPANSION 

We consider from now the massive 0(2N) Gross-Neveu (GN) model [^l], ^2|, though most of the discussion and 
equations are kept general, applying a priori with minor adaptations to other AFT models, as long as the low orders 
RG properties are know. The GN Lagrangian reads 

_ _ o 2 — . 

C GN = -m** + y(#*) 2 (2.1) 

where g 2 is the four-fermion coupling and m an explicit fermion mass. In the D = 2 GN 0(2 jV) model, the discrete 
chiral symmetry \& — ► 75^ of the Lagrangian in the m — ► limit, is spontaneously broken|2l| for any g 2 > 0. 
The model develops a non-trivial mass-gap M P (N), whose exact expression in the massless limit m — > has been 
established for arbitrary N by using the exact S-matrix results [^3| and thermodynamic Bethe Ansatz methods [pi) . 
In this paper we shall mainly (but not only) work in the so-called 't Hooft renormalization scheme where 

p(g 2 ) = dg 2 /d\n^ = -2b a g 4 - 2b l9 4 , (2.2) 

(also defining our RG coefficient conventions, so that bo > for an AFT). This is motivated by the fact that all higher 
order coefficients bi for i > 2 are non universal, being explicitly renormalization scheme (RS) dependent. Similarly 
we truncate the anomalous mass dimension to two-loop order: 

lm{g 2 ) = -d(mm)/dm/i = 7o.g 2 +7i5 4 ■ (2-3) 
In the 0(2N) GN model, the RG coefficients bi and 74 are exactly known in the MS scheme up to three loop order [p8| : 

_ 2V-2 _ 2N-2 ™ (2JV-2)(2JV-7) . 

and 

2JV-1 -ms 2N-1 -ms (2V-l)(4V-3) 
70 ^^^' 71 = ~^^' 72 32^ • (2 ' 5) 

where 70 is universal while, as indicated, the 71 coefficient is RS-dependent, as discussed in more details later (see 
also Appendix ^) . 

A. Pole mass and 1/F expansion properties 

In the above scheme ( |2.2| ),( |2~3| ), one can write the GN pole mass at arbitrary perturbative orders in terms of the 
scale invariant mass rh, as follows [0-10,0: 



M p (m) = 2-° mF- A [C + F]~ B 
with m the scale invariant mass at second RG order: 



1 + J2 d n {2b F)- n 



n=l 



(2.6) 



m = m^)(2b g 2 ^))- A [l + ^g 2 ] B , (2.7) 



F(|) = ln[|] — AlnF — (B — C) \n[C + F], (2. 
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related to the usual perturbative coupling as F 1 = 2&o.9 2 (-^rg) where Mrg is the pure RG resummed mass: 



M RG (m) 



- o-C„ 



F- A [C + F] 



—B 



(2.9) 



A is the basic scale at second RG order in the MS scheme: 



bx 



A = ,ie (6 o5 2 ( Ai ))-c [1+ i^c 



(2.10) 



and the parameters A, B 7 C in Eqs. (2.6)— ( 2.10 ) are given in terms of the one-loop and two-loop RG coefficients: 



7i 
2V 



7o 
26 



(2-11) 



while the coefficients d„ in Eq. (|2.6| ) are essentially made of the non-RG (non-logarithmic), purely perturbative 
contributions from the rt-loop graphs (generically dominant, as will be discussed later), plus eventually (subdominant) 
contributions from higher RG orders in a scheme with 6 i; 7i ^ for i > 2. In the 0(2N) GN model, those perturbative 
coefficients are only known exactly to two-loop order, and indeed similarly in most of the models. One finds[jl4|, ^] 
in the MS scheme 1 : 



' V 2 16 7 2vr 2 



(2.12) 



Eqs. (gg)-(|Dg) were obtained by integrating exactly the usual RG evolution for the (renormalized) Lagrangian 
fermion "current" mass: 



m{fi) = m(/i) exp{- 



9ifl \ 7rn(9) x 



(2.13) 



using the self- consistent condition Mac = m(fi = Mug) defining Mug in Eq- ( |2.9| ). Note that this resummation of 
the pure RG dependence, exac t at second RG order (and to all orders in the two-loop truncated scheme ( |2.2| ), (2.3), is 
explicitly factored out in Eq. ( |2.6| ) from the purely perturbative series d n / (2boF) n : the resummation of the latter 
series is precisely the non-tri vial i ssue in most renormalizable models, as discussed in details later. 



One can easily check that (2.6)-( |2.9| ) are scale invariant expressions, by construction to all orders 2 . Eq. (2.6) is 
perturbatively consistent, for fh A, with the usual expansion relating the current mass at the scale M pote and 
M poie itse if J29I gol: 



M Poie = m (Mf° Ze )[l + J2 c n g 2n (M pole )} 



(2.14) 



where the c„ coefficient are related to the d n ones in Eq. (2.6) in an easily calculable way, involving also RG-dependcnt 
quantities, whose precise expressions are unessential here. 

W e co ncentrate now for the time being on t he p roperties of the pure RG-dependent, resummed mass expression, 
Eq. (Q. At first RG order (61 = 7l = 0) Eq. (|J) takes the simpler form 



F(m/A) = ln(m/A) - A InF = A W[A^ {m/ 'A) 1/A °] 
where the Lambert fl7|l function W[x] = \nx — \nW, is plotted in Fig [ij and 

7o 2N-1 



A = 



2b a 



2N-2 



1 



N 



(2.15) 



(2.16) 



DO 



(where the specific value of Aq for the 0(2N) GN model is indicated for illustration). Eq. (2.15) has the remarkable 



1 We quote in Eq. 
value as given in 



Again 



the exact value of d,2 , recently obtained in the second ref. of 

d,2 ~ (0.737775 — 7r 2 /96) (2./V — 1) / ( 2tt 2 ), does not affect any of our numerical results 
all orders" at this stage means within the scheme fc.2l), (2.3). 



ll . The tiny difference with the former (numerical) 
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FIG. 1: The Lambert W function compared to the Log. 



property: 



F ~ (to/A) 1 ^ (2.17) 

TO — > 



for to — > 0, in contrast with the ordinary logarithm function (see Fig. 1), but the latter is asymptotic to F(m/ A) for 
to 3> A. More precisely, on its principal branch (defined to be the one real- valued for real to), F has an alternative 
power series expansion [[l7[ pp[ : 

which has a finite convergence radius R c = e~ A °(Ao) A ° , of order 0(1): for exam ple in the 0(2N) GN model, R c ~ 0.41 
[i? c = e" 1 - 0.37] for JV = 2 [JV — ► oo]. The pure RG mass M RG (rh) in Eq. ( ^j|) thus exhibits different branches, 
determined both by the values of the RG parameter (— 1) A ° and by the original two branches of the Lambert function 
(see Figs. [TJ and ||). In Fig. || the two principal branches, extending for Re[M]/A < 1, correspond to F < with 
to > or to < respectively, as given by the roots of (— 1) A ° = (— l) 5 / 4 e.g. for N = 3 in the 0(2N) GN model. Each 
of these branches divides in turn into two branches, for e~ Aa A < i?e[Af] < A and < i?e[Af] < e~ A ° A, respectively, 
which correspond to the original branch structure of the Lambert function, above and below W — F = —1, see Fig. [l]. 
Clearly, the physical branch is the one with to > 0, A < Re[M] < +oo. It is indeed the only one branch which for 
real rh values, is real and continuously matching the asymptotic perturbative behaviour of F at large to. Moreover it 
is the branch consi stent with a non-zero " mas s gap" Mrg — A for to — * 0. Algebraically, this RG mass is obtained 



by expanding Eq. ( 2.18 ) to first orders in (2.9) 



M RG (m -> 0) = m [(to/A) 1 ^ + ■ ■ -}- Ao = A (1 + ©(to/A) 1 ^ ) , (2.19) 



which can be viewed as a generalization (for m ^ 0) of dimensional transmutation |18|: more precisely, using Eq. (2. IS) 



one can easily get systematical corrections in powers of (m/A) 1 ^ to the Mrg/A ratio in Eq. ( 2.1 9| ) . Eq. (2.1£) 



automatically reproduces, e.g., the GN O(N) model mass gap in the large N, to — > limit (where A — > 1 for 
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Re[M]/A 




FIG. 2: The different branches of Re[M(m)]/A at first RG order in Eq. ( pjh , determined by the value of A = A in Eq. ( |2.16| ), 
for A 7 = 3 in 0(2N) GN model. Similar branch structures occur at second RG order, for Eq. (2.9) with 61,71 7^ 0, and as well 
in other AFT. 



N — ► 00), traditionally obtained in a different way The multivaluedness of F and Re[M] in the infrared region 
does not signal ambiguities for the physical quantities: actually the different branch structure is completely fixed for 
given RG coefficients. This is generic, so that a similar structure occur s e.g . in QCD, but with of course a different 
value of the RG parameter Aop(J. We again stress, however, that Eq. (2.19) alone is not a proof of dynamical XSB , 
as one may naively infer from the previous discussion: in fact, our construction only exploits that (from dimensional 
transmutation) any mass is proportional to A for rh — > in an AFT, which is encoded here in the properties of 
F(rh) for any rh. In particular we emphasize that the same properties still hold e.g. for any 2-D AFT models with 
a continuous chiral symmetry where spontaneous breakdown is not possible according to Coleman theorem but 
where a physical mass gap occurs pq]. 

At second RG order, F as defined by Eq. (2.8) cannot be written in terms of the Lambert function, but it has very 



similar properties which can be easily inferred from its reciprocal function: 



'j.=e F F A {C + F) B - c . 



(2.20) 



Indeed, replacing Eq. ( [2.20 ) in Eq. ( |2.9| ) immediately gives a simpler expression for the pure RG mass, as a function 
ofF: 

M RG {F) = 2-° e F (C + F)- c A (2.21) 

which thus only depends o n th e universal RG quantity C defined in Eq. ( 2.11 ) . For illustration we plot in Fig || the 
function F defined in Eq. (2.8), for the 0(2N) GN model case with N — 2, w here the relevant RG coefficients bo, 
bi, 70j 7i were defined in Eqs ( [2.2] ), (2.3). Similarly to first RG order in Eq. (2.18), F has a power expansion for 
sufficiently small rh: defining 



n(B-C)IA 
F{x) = AG[- : x 1/A ^ 



A 



(2.22) 



with x — rh/ A, one has from (2 



G = xe- G (l + ^Gf~ B 



P =i 



(2.23) 
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m/A 



FIG. 3: The different branches of F at second RG order in Eq. (2.£) in the 0(2N) GN (N = 2) 



where the expansion coefficients a p are now more involved than the first order corresponding ones in ( 2. lq ) (the a p 
will now depend explicitly on the RG quantities A, B, C), but can be derived systematically 3 . It is easy to determine 
the convergence radius of this expansion form of F(rh) around zero, given by the location of its closest non-zero 
singularity 4 : 



dF y 



and correspondingly 



-> F* 



4 2i 
7o 2 



(2.24) 



!li- = e F' { F*) A {C + F*)( B -^ 



(2.25) 



so that R c = |m*/A|. Note that for N — > oo, F* = — 1 and correspondingly rh* / A = — e 1 . This is due to the fact 
that for N — > oo F is exactly the Lambert function, see Fig. |. At second RG order, in the 0(2N) GN model, F* < 
and C + F* < are real and negative ViV > 3/2, so that the number and location of the singularities in Eq. (2.25) 
on the circle of radius R c are determined by the roots of (— 1) A+B ~ C : see Fig. 0. 



B. Other renormalization schemes 



Before to proceed we sho uld r e mark that at second RG order there is a certain arbitrariness in e.g. Eq. fl2.8|) 
and related quantities, like (2.24)— (2.25), since 71 is renormalization scheme (RS)- dependent: more precisely, for an 
arbitrary perturbative RS change in the Lagrangian mass and coupling parameters: 



g 2 (1 + A l9 2 + ■ 
m (1 + B ig 2 + 



71 is changed as 



71 = 71 + 2b B 1 - 70A1 = 71 + (57i 



(2.26) 



(2.27) 



while bo, b\ and 70 are RS-invariant, as already mentioned (more details on RS changes are given in Appendix C). This 
means for instance that the location of the singularities and the value of the convergence radius as implied by ( 2.25| ) 



3 The a p in Eq . ( |2.23| ) are easily evaluated to high order with e.g. Mathematical?^ . 

4 In Eq. (2.24) the solution 1 — is only valid for 71 ^ 0: if 71 = 0, F* = is not a singularity, as is clear from the first RG order 
result, where the only singularities ly on the circle of radius e~ A ® (Ao) A o . 
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FIG. 4: Convergence radii of the power expansion form Eq. (2.23) of F at second RG order, Eq. (2.8), and corresponding 
singularity locations on the convergence circle for different values of N in the 0(2N) GN model. 



may be mo dified , to some extent, by appropriate changes in 71 (equivalently changes in the quantity A = 71/(261) as 
defined in (fuij)). 

It is always possible to choose the 't Hooft scheme, in which Eqs. (|2.6[) , ( ^.9[ ), ( |2.8|) resum the complete RG dependence 
in rh. This is convenient because beyond second RG order, in an arbitrary scheme where 61,7, ^ for i > 2, algebra 
becomes quite involved, and neither the non-log contributions d n nor the b n and 7„ RG coefficients are known at 
arbitrary orders for most field the ories , and in particular for the GN model. Nevertheless it is still possible to work 
out a formal generalization of Eq. ([2.6|) , in an arbitrary (MS) scheme, see Appendix A. We will use this generalization 



to define some of the numerical approximations to the mass gap, of arbitrary higher orders, in s ectio n VII. 

Before to conclude this section, we discuss another possible RS choice, obtained from expression ( |2.6| ) by an all orders 

redefinition of F: 



F-(A- 



F 

B - C)\n — 
F 



C)M1 + 



(2.28) 



which can be perturbatively expanded in powers of l/F, where F is now again directly related to the Lambert 
function: 



F(^) - Hj] — (A + B — C)lnF = AW[(^ A 



/A] 



(2.29) 



with A = A + B — C. This redefinition is motivated from the fact that in the GN model, the RG coefficient C < in 
Eq.(^TTl), due to bf N < 0, which corresponds to an infrared fixed point at g c = —60/61 > 0, so that the perturbative 
branch of F reaches rh — > first for F = —C. In the scheme (2.28), Eq. (2.6) takes the form: 



M p {m) = (2/e)" c m p- {A+B) ^ d n (26 F)" 



(2.30) 



n=0 



This als o implies appropriate changes in the purely perturbative coefficients, simply determined by re-expanding 
Eq. (O) in l/F powers: 



d\=d\ — B C 



(2.31) 



III. INFRARED RENORMALON PROPERTIES OF THE GN POLE MASS 



As mentioned in introduction, the idea is that, since the complete pole mass Eq. ( |2.6| ) gives the ratio M P (F)/A to 
all perturbative orders for m 3> A, if we are able to resum this series and to give it a meaning for m — > 0, we can obtain 
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the M p /A ratio in the physically interesting massless limit. As far as the pure RG dependence is concerned, this turns 
out to be possible because F{fh) provides a rigorously defined and explicit bridge between the "non-perturbative" 
to < A regime, where F has power expansion ( 2.1§| ), and the short distance perturbative to 3> A (logarithmic) regime. 
A crucial point indeed is the difference between the usual effective coupling g 2 (p 2 ) = l/[ba ln(p 2 /A 2 ) ], having a Landau 
pole at p 2 — A 2 , and i 71- 1 (m) here, having its pole at to = 0, governing the massless limit ( [2.19| ) of the (pure RG) 
mass gap Eq. (2.£). Acc ording ly along the continuous branch on Fig. 2, M(rh) has no singularity for < to < oo, 
as is clear also from Eq. ( 2. IS ) and Fig. 1, 2. Now, to extrapolate the complete pole mass ( |2.6D down to the chiral, 
strongly coupled regime to ~ 0, the main obstacle comes from the presence of the purely perturbative coefficients 
d n . First, though the pole mass (or other physical quantities similarly) is infrar ed fin ite, gauge p9||-, scale- and 
scheme-invariant, the relation between the pol e ma ss and e.g. the running mass in ( 2.14 ) is scheme dependent, which 
is m anifes ted here by the RS-dependence in (2.6) of the perturbative coefficients d n , the RG coefficients A, B in 
Eq. (tZ.llD , and of A too. 

Second, it is immediate that the perturbative contributions d n /F n in Eq. ( |2.6[ ) a rc sin gular when F — ► (m — > 0), 
since each term will have a leading divergence ~ d n (m/ A)~ n / A , according to ( 2.18 ). In other words, while the 
usual Landau pole problem was avoided in the pure RG part ( |2.9| ) of M p , which has a regular finite to — > limit, 
as illustrated in Figs. 1-3, a problem reappears in the perturbative corrections relating the true physical quantities, 
like the pole mass, to their pure RG part. Thus, in an arbitrary scheme, strictly speaking M p — ► oo w hen to — > 0. 
(In principle one could avoid this problem in a crude way by exploiting the RS arbitrariness in (2.6) to define a 
scheme such that all the perturbative coefficients d n = 0. Although such a peculiar scheme can always be formally 
constructed, this solution is to be considered unsatisfactory, since one expects truly non-perturbative results not to 
depend on a particular scheme.) This appears in fact completely similar to the perturbative expansion in powers 
of (g/m 3 ) of the oscillator energy levels, thus also singular for to — ► 0, which nevertheless do not prevent different 
resummation methods to work very well fUjfl , plf , even for to — * 0. We will see in next sections how similar 
resummation properties generalize, to some extent, in the present field theory case. 

Now there is unfortunately an even worse problem, when dealing with the purely perturbative expansion of the pole 



FIG. 5: The GN mass graphs at order l/N. 

mass: in the O(N) GN model, at order l/N, it exhibits infrared renormalons very similar to the QCD quark pole 
mass|Q, . More precisely, let us consider only the naive perturbative expansion of the pole mass, obtained from Fig. 
by taking perturbative expression of the dressed scalar propagator (wavy line), G~ 1 (q 2 ) ~ g 2 {q 2 ) ~ [bo ln(g 2 /A 2 )]~ . 



M 



1 



d^_ 

-LV./o M 2 



(i-0 



M 2 

g 2 { q 2 )-0{\) 
T 



(3.1) 



where £ = (l+4M 2 /g 2 ) 1 / 2 and M = fie V( 2i, o 9 2 M) j s the mass gap at leading l/N order. Then a standard calculation 
givespl 



d„+i „ (2b ) n n 
n — * oo 



(3.2) 



so that the series Eq. ( |2.6| ) including this next-to-leading 1/iV order is badly divergent for any to, and not even Borel 
summable: such a factorial growth of the perturbative coefficients, with no sign alternation, implies Q ambiguities of 
0(A). But, those renormalons are only perturbative artifacts: considering now the full scalar propagator contribution 
(which is known exactly for the O(N) GN model at l/N order): 



G[q 2 ] = 



g 2 N 
2tt 



M 



C + ii 



9 2 N.. r C + L 



Oil/N) , 



(3.3) 
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rather than its truncated perturbative contribution Eq. ( p.l| ), the exact 1/N expression of the pole mass is obtained |24|, 
M as 



M 



A 



In 6 



IE 



21n(cosh[0/2])+ln|^] 



(3.4) 



with X = (1 + 4Af 2 / A i 2 ) 1 /2 = l/tanh(6»/2) (i.e. 6> = ln[( X + l)/(x - 1)] > 0), and £?i(-:c) = -J^dte^/t the 
Exponential Integral function (x > 0). Thus 8] has a factorial perturbative series wit h si gn- alternated coefficients, 
i.e. the IR renormalons actually disappear: more precisely we can re-expand the result (|3.4[) in perturbation, using 



ln[ 



X + l ] 



111 



M 2 



V 



M 2 



9 



(3.5) 



M 



M 



2N 



2 In 2 



IE 



AT 2 °° M 2 

^E(-i)^ 5 2( " +1) +o&)] 



(3.6) 



The explicit Borel summability of the genuine perturbative expansion, Eq.(3.6), is not in contradiction with the 
purely perturbative results above, because the non-trivial cancellation of renormalons involve the contributions of 
non-perturbative power corrections contributions |35|, [3(|. Moreover, it turns out that this cancellation is such 
that the final expression of the pole mass contains neither "purely perturbative" nor "intrinsically non-pe rtur bative" 
contributions: for instance, the net contribution due to the first graph in Fig. [5] is the term —je in Eq. (3.6), which 
simply remains after cancellation of the first order terms: 



dt 



:[1e + 



1 



-t/g 2 \ 



t(l+t) 

M 2 
— ) 



0( 



M 2 



(3.7) 



Now, the point is that the above results Eqs. (3^), (3^) obviously could only be obtained from calculating explicitly 
the exact mass gap at next-to-leading 1/iV order. Our aim here is to ignore on purpose these exact 1/N results, a 
priori only accessible in a certain class of 2-D models. Rather, we want to examine whether our generic construction, 
relying solely on the purely perturbative information, together with the infrared properties of the function F, is able 
to recover some of the non-perturbative properties of the exact mass gap, in particular its Borel summable asymptotic 
expansion explicit at next-to-leading 1/N order. 



IV. BOREL SUMMABILITY OF THE F EXPANSION 



We first reexamine here why the expression (2.6) is plagued with perturbative ambiguities, and how one can get 
rid of those, within our construction, thanks to the analytic properties of F in a vicinity of to = values ]2C[|. First 
we define from Eq. (2J3) its Borel transformed series 5 : 



B.T.[M p {F)][t) = 2- c mF- A (C + F)- B 
so that the corresponding Borel integral reads: 



=j(n-l)l 



BI(m) = M p (m) = 2- c m F- A (C + F)- B / die"* [1 + (4tt b F)' 1 V (^) r 

Jo Z^n F 



(4.1) 



(4.2) 



upon assuming for the perturbative coefficients d n the leading large order behaviour in Eq. (3.2). For any F > 0, this 
expression would be (asymptotically) equal to (2.£) by formal expansion, would the pole at t — F not make the integral 



5 The resummed RG-dependence rhF A (C + F) B , having obviously no factorial behaviour, is thus factored out of the Borel transformed 
series. 
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(4.2) ill-defined. One should make a choice in e.g. deforming the contour above ( or be low) the pole, which results in 
an ambiguity, which is easily seen to be proportional to 0(e~ F ). Since from Eq. ( [2.15 ) F ~ ln[m/A] — A In [In [m/ A]] 
for rh ^> A, this implies a perturbative C(A/m) ambiguity for the "short distance" (M,rh ^> A) pole mass: 



.A Ar m ^ 

ambig ~ ±i ( — ) In — 
in A 



(4.3) 



in consistency with general results^]. Now in our case, Eq. (|2.17|) (and equivalently Eqs . (|2 . 2 2[) , ( |2^) at second RG 
order) allow to trace the behaviour of F all the way down to rh — > + , where at first RG order, F — > 0: consequently 
the naive mass gap (|2.6|), expected to be ~ A, is also ambiguous by O(A). 

But in contrast, within our construction, the Borel integral (4.2) can be defined unambiguously and independently 
of the RS parameter A, in the range F < 0|2(J: then F = —\F\ simply produces the adequate sign alternation in 
the factorially growing coefficients ~ F~ n . More precisely, a straightforward calculation of Eq. (4.2) for Re[F] < 



(neglecting for simplicity at the moment the two-loop RG dependence C, irrelevant to asymptotic properties), gives 



BI{M P /A) ~ , 
F < 



-\F\ 



(4.4) 



where we also used Eq. (2.2C ) to express M p /A as a function of F only. Indeed, as already mentioned the first 
RG order function F(rh) in ( 2.15 ) is well-defined (analytic) for any A values in a disc of radius e~ A A A around 
zero (and for A = 1 the only singularity is at F = —1 i.e. rh /A = -e"\ cf. Fig. 0). Thus, one can choose the 
branch of F such that Re[F(rh)] < 0, compatible with the limit rh — > 0. The second RG order F in Eq. fl2.8|) 
has similar properties, with finite convergence domain around F = and F = —C respectively, see Fig |^ 6 . More 
generally, in an arbitrary AFT with arbitrary values of the RG parameters A, B, C depending on the renormalization 
scheme, there always exist branches of F such that F is complex. This is the case for the GN model for the 
two branches shown with Re[F] < — C in Fig. g|. As a consequence, the Borel singularities in e.g. Eq. (4.2) are 
moved away from the real semi-axis of Borel integration Re[t) > 0, now being located at to = F = \F]e l0F with 6p ^ 0. 



We obtain in this way formal Borel convergence for a certain range of the expansion parameter near the relevant 
massless limit rh = 0, strictly only along those branches such that F < 0, or more generally complex. Depending on 
the branch of F, this may correspond either to rh > or rh < (see Figs. ^|3|), which is in prin ciple not a problem, 
since relativistic field theories only depend on m 2 (the sign of the Lagrangian mass term in ( |2.l| ) can be flipped by a 
discrete 75 transformation, and the Dirac equation is invariant under m — > —m). In Fig. |^, the complex branches of 
F are those corresponding to Re[M]/A < 1, and to i?e[F] < — C in Fig. |3|, where in both cases the symmetry with 
respect to rh — > — rh of such pure RG dependence is manifest. More generally, we expect that our final, physical mass 
gap result, should be independent of the way in which the massless limit is to be reached, either from Re[F] > 0, 
or i?e[-F] < 0, or more generally from any of the complex branches of F. However, the (unphysical) perturbative 
expansion is clearly not invariant under this, since the (usual) expansion with real F > is non Borel summablc 
and ambiguous. Therefore, the following picture emerges: in our construction, the perturbative expansion near the 
strongly coupled, massless limit can exist in two modes: 

i) in the standard mode, corresponding to real F > 0, and matching the usual perturbative expansion for rh 3> A, the 
perturbative expansion alone has to be necessarily completed as usual with "non-perturbative" power corrections, as 
illustrated explicitly with the exact 1/N calculation of the GN mass gap, discussed in section 3. 

ii) In the "alternative" expansion mode, with F < in the simplest case (first RG order) where Aq = 1, the 
perturbative series is directly Borel summable, thus non-ambiguous. There is, therefore, no explicit non-perturbative 
power correction contributions needed in principle. Those results are completelly general for a renormalizable AFT 
of dimension D > 2, since they only depend on the generic RG properties of the function F. 

To illustrate perhaps better the last points, we can give an analogy, to some extent, in the simplest possible model 
where such issues can be discussed, the anharmonic oscillator. The detailed analogy is discussed in Appendix [^. 
The oscillator is described |l3| by a g</> 4 massive scalar field theory in 1-D, with energy levels having from purely 
dimensional considerations a perturbative expansion in powers of g/\m\ 3 : 

Eo~\m\ 5> n (-^r, (4.5) 



Note that at second RG order, from Eq. (2.20) the point F = — C also corresponds to rh = 0: we shall come back on this later on for 
the GN model, where F = —C > corresponds to the infrared fixed point at g 2 = — 6q/6i > 0, as already mentioned. 
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with factorially growing but sign-alternated coefficients |T|, [57|: a n ~ (— l) n n! at large orders, thus the energy levels 
are Borel-summableQ. Now, let us assume, momentarily, that our only knowledge of the oscillator would consist of 
the perturbative expansion Eq. (fh^) , and consider formally changing the sign of the coupling g there: this obviously 
induces a change of sign in the perturbative coefficients, rendering the corresponding series non Borel summable 7 . This 
accordingly produces an ambiguity, an imaginary part in the energy, whose leading terms can be evaluated exactly 
still using the Borel integral, and according to the standard interpretation it calls for additional non-perturbative 
corrections. The latter are easily shown to have the form of the standard instanton contributions to the ground-state 



energy) 37 



ImE 



4 4m 3 .TO 3 , 

e 3 8 

/2?r -g 



™ -1/2 



(4.6) 



where g < 0, and governing accordingly the decay of the wave function due to barrier penetration |37J (see appendix 
[b| for more details). Of course, this instanton contribution was originally not derived from the perturbative ambiguity 
since in the oscillator case a "direct" non perturbative calculation is possible. But as seen from the "perturbative 
only" side, the above argument indicates that, already for the simpler oscillator case, the perturbative expansion can 
exist in two different modes, one in which it is directly Borel summable, while for g < non-trivial non-perturbative 
contributions are needed to get a consistent physical picture. 



Coming back to the D > 2 field theory case, we stress, however, that the above discussed Borel summability 
properties in our construction is possible due to the negative (more generally complex) tail of the perturbative 
expansion parameter 1/F in the infrared, rather than due to an artificial change of sign for any F values. One may 
wonder if such a sign alternation of the badly behaving infrared factorials, may not alter the other way round the 
signs of the UV renormalons, which originally have the good (alternated) signs in AFT[Q. It is easily realized that 
they are in fact unaffected by the infrared properties of F < 0, since by definiti on t he UV renormalons originate only 
from the domain /x ^> A (more precisely the Borel integral equivalent to Eq. (4.2) for UV renormalons corresponds 
to integration from /j, 2 < q 2 < oo). In this range, F is necessarily real positive and larg e, se e Fig. and Fig. § 
Remark finally that, only from the properties of F around F < 0, the Borel sum in (4.4) reproduces qualitatively 
the asymptotic behaviour of the exact 1/N result Eq. (3.4) in the O(N) GN model (with Airbo = 2N — 2), except 
for finite terms je etc, which not surprisingly cannot be guessed by our simple Borel summation of the (leading) 
renormalon asymptotic behaviour in Eq. ( |4.2| ), and with only the first order RG dependence included. Thus, at this 
stage the Borel summability property of the series for F < plays a rather formal role, since the leading order Borel 
sum Eq. (4.4) is not expected to be a good numerical approximation of the exact mass gap. (We shall see in section 



VII that there are more efficient approximations to the exact mass gap.) 



It is not difficult to work out the exact second RG order generalization of Eq. (4.4). For this purpose, it is convenient 
to define another change of scheme, by a "non-perturbative" redefinition of F: 



F + C = f 



(4.7) 



perturbatively equivalent to 



5 2 



g 2 (i + jvr 1 



(4.8) 



This redefinition is again motivated from the fact t hat in the GN model, the RG coefficient C < in Eq.(2.11), so that 
the extrapolation from the perturbative range of (2.9) down to to — > reaches first F = — C, as explained before(see 
Fig. ||). The asymptotic behaviour of the perturbative coefficient (exact in the scheme ([T^)) is 



dn+l 

71 — > OO 

for n — > oo. The corresponding Borel integral reads 

BI{M p )[f]=K{2e)- c eff- C [1 



(2b ) n T[n + l + C] 



r[i + C] 

47TD / 



die'* (1 - t/f) 



-l-C ■ 



(4.9) 



(4.10) 



7 This is not to be confused with the double- well potential, obtained from the oscillator by the change: m 2 — > — m 2 , which also has a non 
Borel-summablc perturbation series. 
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where we used the appropriate RS change Eq. (4.7) in Eq (2.6). So for F < —C, i.e. / < 0, we obtain 
BI(M p )[-\f\]~A(2e) 



- C e-^Hf\r C [1 - A" T[l + C] e\f\ \ff T[-C, |/|] 



(4.11) 



where the 1 in the bracket refers to the pure RG part, while the remaining part in Eq. ( 4.11 ) resums the purely 
perturbative contributions. The latter resummation of the perturbative expansion part is convergent, giving a finite 
contribution V|/| ^ oo. In particular it gives a finite result even for |/| — > 0: 



BI[pert. series] 



1/1-0 



(~l)- c (2eY 



47T b 



r[i + c] r[-c] 



(4.12) 



which accordingly, corresponds to reac h the massle ss limit m — > along any of the branch of F below F = —C, see 
Fig. |[ Remark also finally that Eqs. (4.11), (4.12) only depend on the universal RG coefficient C in (2.11) (noting 
also that 47r& = -C _1 = 2N - 2 in the GN model). 



V. VARIATION ALLY IMPROVED MASS EXPANSION 



We examine now how to complement the above construction, based only on RG fixed point properties, by combining 
it with a specific variant of the delta-expansion, or variationally improved perturbation (DE-VIP) method. The latter is 
us ually[f§ |l§ applied more directly on the ordinary perturbative series in the coupling g. But the present construction 
and the DE-VIP idea are closely related, the DE-VIP being essentially a reorganization of the interaction terms of 
the Lagrangian, with the introduction of a trial mass parameter, for physical quantities relevant in the massless limit 
of the theory. We will see that the DE-VIP can give further improved Borel convergence properties. We define the 
(linear) DE as the power series obtained formally after the substitution 



(5.1) 



within any perturbative expressions at arbitrary order, where to(ju) is the renormaliz ed L agrangian mass (in e.g. 
MS scheme), 8 the new expansion parameter, and m v an arbitrary adjustable mass. ( |5.l| ) is equivalent to adding 
and subtracting to the massless Lagrangian a "trial" mass term m v [8 interpolating between the free (8 — 0) and the 
interacting massless Lagrangian (8 = 1)], and is entirely compatible with renormalizationjlil and gauge- invariance ]l6[ . 
After substitution (5.1), 



M p (m,8) = ^a k {rh)8 k 

k 



(5.2) 



can be most conveniently directly resummed, for 8 — > 1, by contour integration |14j around 5 = 0, to ar bitra ry order 
K: an appropriate change of variable allowing to study the m(/i) — > (equivalently 6 — > 1) limit in Eq. (5.1) is: 



8= 1-v/K; 



= K 1 



(5.3) 



Eq. (5.3) is simply a convenient way of parameterizing how rapidly the Lagrangian mass m(fi) — > limit is reached 
(as controlled by 7 < 1) as function of the (maximal) delta-expansion order K. Similarly to refs. O, 12 the point 
is to adjust the rates at which m(/i) — * (8 — * 1) and K — > 00 are simultaneously reached, with no a priori need of 
invoking explicit PMS optimization principle. The final contour integral summation takes a simple form, for K — > 00: 



M p 
A 



,-c 



1 

2~7ri 



d ve Wm") F- A [v] (C + F[v])- B [l + 



1 



N 



2b F[v}^ o F^[v\ l 



(5.4) 



where m" = m„/A, N is the maximal perturbative order, and after a deformation t he c ontour encircles the semi-axis 
Re[v] < (see Fig. |^). For simplicity we fix from now the scaling parameter in E q. (|5.3| ) to its maximal value (7 = 1) 
still compatible with massless limit (i.e. m v — * 0). (NB The general 7 scaling (5.3) can be analyzed in a way more 
similar to the oscillator |ll], ^2| , i.e. without the peculiar contour (^-summation Eq. (5.4), see Appendix D. But this 
gives rather cumbersome algebraic and numerical analysis.) Eq. ( [5.4| ) can be well approximated analytically (at least 
for slightly restricted RS choices, as will be indicated): 



M P /A - 1 + 



47T&0 



N 

E 

9 =1 



N-q 

E 

p=0 



T[p + q}(p + q + A)(q + A)P- 
APT[l+p}T[l + q/A} 



(m") 



-q/A 



(5.5) 
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FIG. 6: Singularities and equivalent integration contours in the v plane, for A = 1. 



where we used Eq. (2.18), leading to the (exact) expression 

F 

together with 



-n-A, \ ^(n + A){n + A-p)P-\ u 



p=0 



p\ Ap 



~~~~~7 I dve v v z = —r- = 
27r»J T[-z] 



(5.6) 



(5.7) 



and also assumed the leading renormalon behaviour 8 Eq. (ET2). In Eq. (5J5) we also made a convenient reshuffl ing o f 
summation indices, n — > p + q, where n is the original perturbation order, p is the order of the expansion in Eq. (2.18), 
and q is the order of the (resulting) expansion in (m")~ q / A . 

In fact, some restrictions apply pCj to Eq. (|5.5|): the sum over p is bounded as given, iff 



i/A e N* 



(5.8) 



since l/r[l + (n—p)/A] = for p > n + 1, which we assume in this section for simpl icity . (This is not much r estric tive, 
except that for arbitrary AFT it is generally not possible both that A s atisfi es (5.8) and B — in Eq. (2.11), as 
assumed in ( |5.4|) . But the more general scheme B ^ simply makes Eq. (5.5) algebraically more involved, without 
affecting the larg e or der behaviour and convergence properties.) 

Second, strictly (5.5) is valid only asymptotically, for sufficiently large N: due to the finite convergence radius of 
expa nsion (2.18), interchanging the sum in ( 2.18| ) and integration in (5.4) is not rigorously justified. However, when 
(5.8) holds, the formerly branch point v = is simply a pole, which allows to choose an equivalent contour of arbitrarily 
small radius aroun d v = 0, thus always inside the convergence radius of (2.18) (see the dashed small circle contour in 
Fig. 3). Using Eq. (5.6) (exact for \m v v\ < e~ A A A ); and e v = J £2 k v k /kl, it is easily seen that the simple poles are for 



k - (1 + (n - p)/A) 



-1 



(5.9) 



for which the contribution to the contour integral will be the coefficient ( 5.6 ) , divided by k\ : this give r [l + (n—p) /A] = 
r[l + q/A]. So, only the simple pole terms i> -1 contribute to Eq. (5.4), which sum up to give Eq. (5.5) again. The 
extra contributi on (around the cut at v = — e _1 , e.g. for A = 1) gives the difference between the "exact" integral (5.4) 
and expansion ( |5.5| ), and can be evaluated numerically. These contributions are easily shown for A — 1 to contribute 
as 0(e~( em ' )h[N] relative to (5.5), where h[N] rapidly decreases for N — > oo. In Table [j] we compare, for the first 
20 values of the perturbative order n, and for m" = 1, the exact integrals: 



jexact = _ I dve v/m" '^-n-Ar 1 

2ni J 11 



(5.10) 



The original nl coefficients in Eq. (3.2) correspond to T[p + <; 
easily be implemented: it essentially replaces (n — 11 ! — > F[n 
the convergence properties discussed below, cf. Eq. (4.9). 



in (5.5). Higher order refinements on infrared renor 
- C] (1 + ri (scheme) / n + ..) where n depends on RS 



alon structure may 
without affecting 
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n 


I e n xact {m" = 1) 


7* e " es (m" = 1) 


-e 1 cut contribution 





1.137545624095354 


1. 


0.1375456240953545 


2 


4.992495983527079 


5. 


-0.00750401647292076 


4 


8.165946096293188 


8.166666666666666 


-0.0007205703734776136 


6 


8.780755015889008 


8.780555555555555 


0.0001994603334534161 


8 


7.28191740833204 


7.281944444444444 


-0.00002703611240395531 


10 


5.01821660133935 


5.0182137345679 


2.86677144956 10 -6 


12 


2.99921746894628 


2.999217731214259 


-2.62267978889 10" 7 


14 


1.597938546376994 


1.597938525158887 


2.121810749323 10" 8 


16 


0.773538637086908 


0.773538638583302 


-1.4963940087 10" 9 


18 


0.3450035964997671 


0.3450035964141139 


8.5653262261 10~ n 


20 


0.1432800076663147 


0.1432800076691064 


-2.79173906215 10~ 12 



TABLE I: Comparison between Eq. (5.10) (evaluated numerically), Eq. (5.11), and the extra cut contributions 



which we calculated numerically with Mathematical 32 , with the approximation using the expansion ( |5.6| ) (thus 
neglecting the extra contribution around — e _1 ): 



A (n + A)(n + A-p)P~ 1 
2^p\APT[l + (n-p)/A] 



(m v ) 



-l+(p-n)/A (5-11) 



the extra e 1 cut contribution, evaluated 



that leads to Eq. ([T|). We also give in the third column of Table 
independently numerically for consistency. 

When m" — ► 0, the discrepancy between the exact integral and the analytical resummation in Table | decreases 
rapidly, as expected, even for small n — > (but the numerical integration becomes unstable for very small m"). If 
A =/= I and B arbitrary, contributions from extra cuts are not so simply estimated, and we were only able to check 
numerically that they are negligible with respect to (|j|) for sufficiently large N.] Thus for large enough N (and/or 
small m") those contributions are unessential for the convergence properties discussed below. 

The factorial damping of coefficients, as compared to the original perturbative expansion, is explicit in Eq. ( [3.5] ). 
Yet, the damping is insufficient to make this series for N — * oo readily convergent. For any low p <C N, renormalon 
factorials are overcompensated if A < 1, but the T[l + q/A] damping decreases in str engt h as p increases, giving 
increasing contributions to the sum over p. The leading contributions to the coefficients of (5.5) happen at intermediate 
values of p. Nevertheless, t he id ea of damping factorials from appropriate RS choice does survive, when considering 
the Borel transform of Eq. (|5.5|), as examined next. 



VI. BOREL CONVERGENCE OF DE-VIP 

We are now ready to combine the previous DE-VIP behaviour of the series in to/ A with the general Borel convergence 
properties for Re[F] < of the initial perturbative series in 1/F, which were examined in section 4. For completeness 
we consider both a linear and non-linear version of the DE-VIP construction, also for reasons that will be clear below. 
We will see that the DE-VIP expansion can generally improve (accelerate) the Borel convergence properties of the 
1/F expansion. 

A. Linear method 



For any given choice of contour avoiding the pole in the Borel plane t (or cut at higher RG order, see Eq. ( [1.1 
can apply the DE-VIP as defined i n sec tion 4, in trodu cing the ^-expansion and contour resummation as in Eq. (| 
but now on the Borel integral Eq. (4.2) (or Eq. (4.10)). It leads to: 



dve v 



2iir J F A (C + F) 1 



die' 1 [1 



1 



4tt b F 



(l-t/FY 



(6.1) 



16 



where F = F[rh v] and the integrand is to be be understood as its formal expansion in t/F[rhv]. Inte rcha nging the 
contour and Borel integrals, the Borel transform integrand, which is a function of i, is essentially Eq. ( |5.5| ) but with 
the replacement 



(6.2) 



standard from the B orel transform, except that in our case we did some reshuffling of summation indices, as already 
indicated after Eq. J5.5| ). One can find after some algebra the asymptotic behaviour for (the original maximal 
perturbative order) N — > oo: 



M^ ar (m") ~ const. A [1 + 



' dt 

47T b 



E 



T[l + q/A] 



(6.3) 



wh ere w e work again here for simplicity at first RG order, neglecting the pure RG res umm ed (C + F)~ B term 
in (4.2). It thus appears that the asymptotic behaviour of the Borel integrand in Eq. (6.3) is that of an entire 
series (at least for A > 0), i.e. with no poles for < t < oo. More precisely, the pole at to = 1 in the original 
(standard) Borel integrand has been pushe d to to ^ +°o due to the factorial damping, so that the Borel integral 
is no longer ambiguous. However, integral ( |6.3[ ) is badly divergent at t — > oo, at least for Re[m"] > 0, so that the 
series is not Borel summable for standard (perturbative) m" values. In fact it is i mpor tant to remark that this 
non convergent result is obtained when considering the exact expansion of F, Eq. ( 2.1S ), to all orders. Naively, 
one may have thought th at th e contour integral would be dominated for q ~ N 3> p, by the apparently " leading" 



terms for m"v c± in Eq.(5.4). But the asymptotic behaviour of our complete series Eq^JSJ)), as well as its Borel 
transform Eq. (3.3), appear much less intu itive than e.g. the oscillator energy levels expansions. Therefore , it is 
due to the "late" terms of expansion ( 2.18 ) (corresponding to the terms p ~ N — q when N — > oo in Eq. ( |5.5|) ), 
that the Borel integral ultimately diverges. We have checked by direct numerical contour inte gration of the 
Lambert function, which can be done e.g. with Ma them atical! , that as N — > oo expression ( |5.5| ) is asymptoti- 
cally correct (see Table |), while a finite truncation of (2. IS) to the first few terms would not be a good approximation. 



Now conversely, the integral in Eq. ( |6.3| ) can converge, for Re[m"] < 0. This is th e case at least for A = 1, which 
can always be chosen by an appropriate and simple RS change, according to Eqs ( 2.27 ). (In particular, such RS 
change only affect the very first ord inary perturbative order coefficients (see Appendix C), thus cannot modify their 
behaviour at large orders, Eq. (3.2)). Now, since m" = m v /A is an arbitrary parameter, it should be legitimate 
to reach the chiral limit m" — ► 0, of main interest here, within the Borel-convergent half-plane Re[m"] < 0. For 
A =^ 1, one may also choose the arbitrary parameter m" with i?e[(m") 1 / j4 ] < such that (3.2) converges, though 
this appears not always possible for any arbitrary A values. Yet, the general Borel convergence properties obtained 
in section 4 prior to the DE-VIP transformation of perturbative ex pans ion, was valid independently of A values. In 
fact, the scheme choice limitation in obtaining convergence of Eq. (6.3) lis only an artifact of our simplest choice of 
the 5-expansion summation defining the DE-VIP series and leading to (|6.3|), as examined in next sub-section. 



B. Non-linear variational expansion 



We consider now a convenient modification of our variational expansion method, defining a non-linear transformation 

First we remind that Eq. (2.15): 



in the variational parameter m„, as compared to Eq. (5.1 



F = ln^ 
A 



AlnF 



(6.4) 



defines F as a systematic power expansion in (m l ,/A) 1 /' 4 , cf. Eq. ( p. IS] ), rather than a single power of m v /A, as 
reminiscent of the renormalization logarithms. This is the main difference (and main source of complexity) with the 
simpler oscillator energy levels, for which each perturbative expansion order is a simpl e (g /m 3 ) n power fll3[. Now, 
since we introduce a reparameterization of the Lagrangian interaction terms, cf. Eq. ( |5.l| ), it may be possible to 



9 Als n hv si rnpl n analogy with the oscillator, which has a series in g/m s . 
Eq.( |2.lj ),(5.5), and the replacements A -> 2/3, m" -> m 2 g -2 / 3 . 



cf. Eq. (Q), thus formally identical to taking 



in 
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remove the logarithm from F by an appropriate definition of the alternative interaction terms. Heuristically, this can 
be achieved by redefinition of the arbitrary mass parameter, for example: 

rh v — > rh v e m v = rh' v (6-5) 
which immediately gives, replacing in Eq. ( |6.4| ): 

F = (m' v y/ A (6.6) 
i.e. a single power dependence on rh' v . Of course, the above manipulation i s jus t an equivalent change of variable: the 



theory is completely equivalent in terms of m' v , since all the complexity of (6.4) is now hidden in the relation between 



m v and to', . Nevertheless, if one can make suc h a transformation to depend on the DE-VIP expansion parameter 5 



in Eq. (5.1), the contour integrand in Eq. ( |5.4| ) resummi ng t he DE-VIP expansion can have a simpler ^-dependence. 
More precisely, instead of the linear DE defined by Eq. ( |5.l| ), we consider in the (variationally modified) Lagrangian 
a mass term 

(l-5)m v g(m v ,5)Vy (6.7) 

where Q now explicitly depend on the new DE-VIP expansion parameter 5, thus the non-linearity. Q{m v ,8) is only 
required to be a scale (RG) invariant function (so that it does not affect any of the RG properties) and to have a S 
dependence only constrained by requiring that it still interpolates between the free massive Lagrangian (for 6 — > 0) 
and the interacting massless theory (for 6 — 1), the original theory. The function Q thus defines a (non- linear) 
interpolation function which is otherwise essentially arbitrary. This may be viewed as a particular "order dependent 
mapping" (ODM)§. 

Specifically we will use here a convenient 10 form for Q: 

m v /A = to" -> m'g(m",v) (6.8) 

G = e F F A -\C + F) B -° (6.9) 
with F = F[m" v], i.e. a deformation of the mass term depending also on the RG parameters. This gives 

F = ln(m"v) — AlnF — (B — C) ln(C + F) = ln(m'u) + F - InF (6.10) 

O F = (m'v) (6.11) 

With the simple power form of F in Eq. ( |6.1l[ ), the contour integral Eq. (5.4) now reads: 



M p _ n 1 fdv ' 



2-c 1 I ^ M( C + to' v)- c [1 + —L— ^ ] (6-12) 



A 27T i / v 47T 6n to' v (to' v) r 

J u n=0 y ' 



which is immediately integrable (again neglecting here the higher RG order C dependence for simplicity): 

n— 



MP - 1 | 1 V ( m ' + 1 ) n+l f613) 



which is now valid for arbitrary values of the RS parameter A. With d n ~ n! for n — > oo the series ( pla| ) thus 
converges iff 

I^Kl (6-14) 



The choice of the non-linear interpolation is clearly not unique. Other choices J3!j] differ in the explicit form of the resulting DE-VIP 
expansion, but have similar asymptotic properties at large perturbative orders. 
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which is only possible if — oo < Re[m'] < —1/2, thus for negative values of the arbitrary mass m', in consistency with 
the linear method results in previous sections. Taking more specifically the asymptotic behaviour of the perturbative 
coefficients Eq. (|3.2|) at first RG order, the series can be resummed to 



M p 
A 



1 



1 



47T b 



ln[l 



1, 



(6.15) 



In summary, we obtain from this non-linear method a directly convergent series, where the convergence domain is in 
the Re[m'] < ra nge o f the trial mass parameter, and this result is now independent of the A values, as expected. 



The series in Eq. ( |6 . 13 ) may also be straightforwardly Borel resummed: this gives 



M 1 



Afl- 



47T 6n 



dt 



(e 



-i(l-r) 



(6.16) 



with r = (to/ + l)/m' , which converges to Eq. (3.15) for Re [r ] < 1, thus for Re[m'] < 0. Since the series in Eq. ( 6.13| ) 
was already convergent in the m! domain given by Eq. ( |6. 14 ) , the role of the Borel summation in this case is simply 
that it can extend the convergence domain, namely from (6.14) to the full half-plane Re[m'] < 0. Another advantage 
of such non-linear transformations is that the extra singularities, at v = e~ A (— A) A , are absent in this picture (in fact, 
they have not completel y di sappeared, no more than the RS A-dependence, but all this is hidden in the relationship 
between to and to' Eq. (6.8), by definition). The numerical application of this non-linear method will be illustrated 



in the end of the next section. 



VII. NUMERICAL RESULTS 



We are now in a position to analyze and compare in some details the numerical results obtained from different 
possible approximations to the 0(2N) GN model mass gap, defined from our construction. The Borel (or direct) 
convergence properties as obtained in sections 5 and 6 are encouraging but unfortunately only formal pr ope rties of 
the large orders, asymptotic behaviour of the true series. For instance, the leading order Borel sum Eq. (4.4) is not 
expected to be a good numerical approximation of the exact mass gap, in analogy with the Borel summable oscillator 
case, where the direct Borel sum expression (see e.g. Eq. (B4) in Appendix B), is of not much practical use for an 
accurate numerical determination of the energy levels. (Apart from a direct numerical solution of the Schrodinger 
equation, the most efficient analytical methods to determine the oscillator energy levels accurately are either Pade 
approximant techniques^], or the ODM [|| or related optimized delta-expansion |ll|, [l2| methods). In D > 2 field 
theories it is even less expected that the first few perturbative coefficients (recalling that only the first two are exactly 
known in the GN model) are close to the asymptotic behaviour. The precise numerical values of the first perturbative 
orders could thus be a priori of much relevance for a precise determination of the mass gap, in particular for low N 
values. 



A. Direct mass optimization 



We start in this first sub-section by exploring a numerical approximation more directly motivated by the usual 
DE-VIP or related "principle of minimal sensitivity" (PMS)|10 ideas. Since the main parameter in our construction 
is the arbitrary mass rh v , and the exact result is in the massless Lagrangian limit, thus independent of the mass, 
the PMS leads naturally to optimize our expressions for the mass gap, e.g. (5.4), with respect to this trial mass 
parameter. We can compare these optimization results when taking successive orders of the original perturbative 
expansion into account. Numerical optimization results are summarized in Tables O a nd III. More precisely, in Table 
H we look for extrema in to/ A = to" of the original mass gap expression Eq. ( |2~6[ ), truncated to first and second 
perturbative orders, respectively, in the MS scheme. In other words, only the exactly know perturbative information 
is taken into account in these results. (Note in particular that none of the above discussed large order, resummation, 
and eventual Borel convergence properties are taken into account in any way he re.) I n prac tice we rather use the 



change of renormalization scheme as explained in section 2, which leads to Eqs. (2.30)— (2.31), more appropriate to 
the GN model. As one can see, there exist optima in to" for all cases studied, but the optimal M p /A mass gap 
results are rather far away from the exact ones, except when N — > oo (where the correct result M p /A = 1 is always 
recovered). Moreover, when adding more perturbative coefficients, there is even a substantial degradation in these 
optimal results. This may appear a bit surprising, if comparing with the excellent results obtained from a similar 
optimization with respect to the Lagrangian mass of the oscillator energy levels Q. In our opinion it simply reflects 
that in the more complicated field theory case, a naive application of PMS ideas may not always work so well, in view 
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of the numerous problems that afflict the original perturbative series, as emphasized in previous sections. 
Next, we show in Table [II the results from a similar optimization with respect to m" , but this time on the DE- 



N (0(2N) model) 


exact 


order 1 (pure RG) 


order 2 


2 


1.8604 


2.215 


3.28 


3 


1.4819 


1.730 


2.63 


5 


1.2367 


1.464 


2.11 


8 


1.133 


1.337 


1.84 


oo 


1 


~ 1 


~ 1 



TABLE II: M p /A from direct optimiz ation of Eq. (|2 
with respect to m", in the scheme Eq. ( 



(truncated to first few orders, with exact perturbative coefficients) 



VIP resummed, contour integral expression of the mass gap: this is essentially Eq. fl5.4D , except that again only the 
exactly known perturbative coefficients are taken into account, thus truncating Eq. ( |5.4|) at first and/or second order 
respectively (as indicated in Table III). The results are much better than those in Table |]|, which may be attributed to 
the expected improved properties of the contour integration resummation. In particular, the results at first order are 
very close to the exact mass gap, at least for N > 3. We observed that for low N values the optima is for rather small 
0.087 for N — 2 in the second column of Table III). When N increases, m" t also increases, up to 



"■opt 



(e.g. 



b opt 



a maximum m'^ pt ~ 0.26 for N ~ 6, and then again vn!' opt — ► as N — > oo. In connection with our Borel convergence 



results of section 4 it is interesting to remark that there also exist optima for m" < values, see Table IV. However, 
inclusion of higher perturbative or ders makes the results to degrade rather rapidly, contrary to what could have been 
expected, as one can see in Table III. We conclude that such a rather naive optimization, keeping only the lowest 
perturbative orders into account in direct inspiration of the PMS/DE-VIP ideas, is not much conclusive within our 
framework. Note also the excellent optimization resul ts obtained, if considering only the first order, i.e. the pure 
RG dependence only. (In fact, in the scheme Eq. ( 2.3C ), the "first order" in Tables UJ, III onl y in volve the pure RG 
dependence, since in the O(N) GN modeljlj] the first perturbative order coefficient in Eq. ( |2.6| ): di (MS) = 0, so 
that d% 



-BC from Eq. (2.31) 



N (0(2N) model) 


exact 


order 1 (pure RG) 


order 2 


2 


1.8604 


2.0493 


2.85 


3 


1.4819 


1.5156 


2.10 


5 


1.2367 


1.257 


1.62 


8 


1.133 


1.152 


1.40 


oo 


1 


~ 1 


~ 1 



TABLE III: M p /A from direct optimiz ation of Eq. (5.4) (truncated to first few orders, with exact perturbative coefficients) 
with respect to m", in the scheme Eq. (b 



N (0(2N) model) 


2 


3 


5 


8 


oo 


M opt /A 


1.929 


1.357 


1.106 


1.024 


~ 1 



TABLE IV: M F /A similarly to Table III optimizing with respect to m", where m" pt < 0. Only the first order perturbative 
contributions are shown here. 



B. Pade approximants including 1/N 2 contributions 



The second kind of approximation that we have performed, is still not related to the Borel convergence properties as 
described in the above sections , be in g ag ain essentially based on the original basic perturbative expansion, but treating 
the variational expansion in (5.1 )— (5.4 ) in an alternative manner to extrapolate in the infrared, non-perturbative 
region. More precisely, in order to define the massless limit m" — > without resumming the complete series, we have 
used a variation flif of Pade approximant || (PA) techniques. The latter are know to give a reliable resummation 
procedure of perturbative expansion in various situations, provided that the singularities of the original expansion 
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are controllable to some extent. The detailed construction of our particular PA is given in Appendix || These 
approximants have the generic form 



const. Aq exp 



1 



du P[ Pj9 ] (u) 



du i P [p,q] ( U ) ~ 1 



C + l. 



(7.1) 



where, prior to the PA con struction, the very same contour integral as the one in Eq. (5.4) has been performed (see 
Appendix ^J). In Eq. (7.1) the overall constant takes into account any terms from e.g. RG dependence (such as 
typically the 2~ c factor in Eq. ( |5.4| ), or depending eventually on the choice of renormalization scheme). The PA 
functions P[ Pig ](u) are, as usual||, rational fractions of polynomials in the relevant variable 1/u ~ 2b g^tt, which is 
related to m" in Eq. (5.4) as follows: 

\B-C _ m ii 



u A e l 



(C- 



(7.2) 



and accordingly have by construction the property of disentangling the usual perturbative In [ln..[m"]] behaviour. 
Standard perturbation th eory corresponds to u — > oo and the massless limit correpsonds to u — > 0. The splitting of 
integration ranges in Eq. (7.1) is due to a necessary subtractio n of UV divergences, uniquely fixed by the perturbative 
expansion. By construction (see Appendix ^J) the PA in Eq. (7.1) have a finite, regular u — > limit, which somehow 
selects a limited (but not unique) form of possible PA, for a given order of the original perturbation. 
In ref. [M, using for the 0(2N) GN mass gap only the first two orders of the original perturbative series coefficients 
in Eq. (EJj), we had obtained from appropriate p, q orders of the Pade approximants, reasonably good numerical 
approximations of the exact mass gap|M: 



Ml 



A 



= (4e; 



l/(2JV-2) 



/r[l-l/(2JV-2)] . 



(7.3) 



The advantage of these PA is that we can systematically incorporate higher orders (approximations) of the original 
perturbative series. Accordingly, to check the stability, and eventual numerical convergence in our approach, we 
incorporate in this analysis a definite information on higher perturbative orders (recalling that their exact expressions 
beyond two loops are still unknown in the GN model). To this aim we exploit the exactly known 1/N and 1/N 2 
dependence of the GN model RG coefficients, obtained in the MS scheme from an analysis near the perturbative 
critical point in a 2 — D = e-expansion in ref. |3q| . The idea is that, since all the perturbative coefficients are RS 
dependent (as discussed in section 2 and Appendix [d]) one can transfer the (exact) RG information of order 1/N 2 
into perturbative coefficients d n of order n > 3 simply by an appropriate RS change. This RS change is essentially a 
change from the MS scheme, in which the e xact 1 /N and 1/N 2 RG dependence was obtained [Q, to the "two-loop 
truncated" 't Hooft scheme, defined in Eqs. (2.2), ( |2.3|) . More precisely, beyond second order in an arbitrary scheme 
where 6j, 7$ 7^ for i > 2 (such as the MS scheme), the pure RG dependence can always be expanded in perturbation, 
thus taking the form of specific contributions to the coefficients of 1/F n in Eq. (2.6) (see Appendix A for details). 
Actually, this does not generate the exact 1/N 2 perturbative coefficients, but only an approximation of these, in a 
certain scheme. In order to uniquely fix those resulting perturbative contributions, coefficients of 1/F n , we also need 
to fix the perturbative coefficients in the original (truncated) scheme, that we assume for simplicity to be zero. This 
may be viewed again as a particular scheme choice. Whether this assumption is a good approximation or not can only 
be decided by the numerical analysis. The precise link with the 1/N and 1/N 2 information, and how the RS change 
generates perturbative coefficients is detailed in Appendix After some straightforward algebra, we obt ain in this 
way systematic corrections to the mass gap in the form of arbitrary order perturbative coefficients in Eq. ( |2.6[ ). This 
gives e.g. for the first few order terms: 



di 



(N-l/2) 
(N — l) 2 
(N- 



-0.3025 



1/2) (iV - 1.4708)(iV - 0.0127) 



0.9375 



(N - l) 4 
(N — 1/2) (AT + 0.1645)(A^ - 



0.6986)(A^ — 1.1486) (AT — 1.3595) 



(N — l) e 



(7.4) 



etc. Including these results into the PA Eq. ( 7.1 



we obtain the numerical results shown up to fourth perturbative 
order in Tables [v] and VI for two different RS choices below. Since by construction those PA are defined in the 
massless limit, the only remaining arbitrariness is the one due to the scale (or scheme) dependence 11 . This leads us 



11 Defining the RS change from the MS to the truncated scheme is uniquely fixed except for one RS parameter, see Appendix 
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to study also this remaining scale dependence by looking for possible optima|TJ] with respect to an arbitrary scale 
change: 

H —y a /i . (7-5) 
From Tables [v] and [vj we can see that there are indeed always optima in a, and we observed that, as the order 



N (0(2N) model) 


exact 


pure RG 


order 2 


order 3 


order 4 


2 


1.8604 


2.1213 


1.6206* [1.524] 


2.044 [2.0434] 


2.1103 [2.097] 


3 


1.4819 


1.4865 


1.3456* [1.344] 


1.5238 [1.4914] 


1.5372 [1.4986] 


5 


1.2367 


1.2268 


1.1917* [1.186] 


1.27375 [1.2397] 


1.2733 [1.2394] 


8 


1.133 


1.1258 


1.1205* [1.109] 


1.16518 [1.1355] 


1.1631 [1.1344] 


oo 


1 


1 


1.0165[0.9972] 


1.0165(0.9972] 


1.0165[0.9972] 



TABLE V: M F '/A from a Pade approximant P[2, 3], in the scheme defined by Eqs (2.3C )- 



respect to the scale parameter a. The order 2 values indicated as *, were obtained in ref. ]14 
values (i.e. a = 1) are indicated in brackets. 



(ffT^) , (ff75|) , optimized with 
Also, the (unoptimized) MS 



increases, these optima become relatively flat, and are closer to the original MS scheme (corresponding to a = 1). 
Those properties are empirically quite satisfactory, as one expects on general grounds that the approximations should 
be less and less sensitive to the scale (or scheme), since the exact result does not depend on the latter. Moreover, the 
comparison of the second, third, and fourth order indicates a reasonable numerical agreement with the exact results. 
Indeed the third and fourth orders in the MS scheme [the numbers in brackets] are very close to the exact results, at 
least for N > 3. 



We also give in Table VII, for indication, the numerical values of the perturbative coefficients up to fourth order 



N (0(2N) model) 


exact 


pure RG 


order 2 


order 3 


order 4 


2 


1.8604 


2.1213 


1.872* [1.707] 


2.079 [2.057] 


2.048 [2.008] 


3 


1.4819 


1.4865 


1.487* [1.486] 


1.528 [1.519] 


1.510 [1.504] 


5 


1.2367 


1.2268 


1.265* [1.252] 


1.270 [1.251] 


1.262 [1.245] 


8 


1.133 


1.1258 


1.163* [1.145] 


1.162 [1.140] 


1.157 [1.137] 


oo 


1 


1 


1.0165[0.9972] 


1.0165(0.9972] 


1.0165(0.9972] 



TABLE VI: M p /A from a Pade approximant P[2,3], in the original scheme corresponding to Eq. (|2 
captions as for table |v| 



Otherwise same 



that enter the PA analysis, for both schemes corresponding to the results in Table ^ [and |V]| , respectively. One can 
see for instance that the perturbative coefficients for N = 2 are increasing as the perturbative order increases, while 
it is the reverse for N > 3, which may explain why the results become better for N > 3. 

A more curious and perhaps rather remarkable fact, is that the PA results are also very close to the exact mass 
gap, w hen considering only the pure RG dependence: i.e. neglecting all perturbative orders, in which case the PA 
in Eq. (7.1), after the appropriate subtraction of divergences (see Appendix |c|) , essentially reduces to a constant de- 
pending only on the RG parameters A, B, C. This is indicated as "pure RG" in the third column of Tables and [v| 
If taking those PA results at face value, it seems that the pure RG approximation gives very good results, for N > 3, 
then there is some degradation when only the lowest perturbative orders are included, and then results become again 
closer to the exact ones when more perturbative orders are included. At this level, we cannot completely exclude that 
this behaviour maybe a numerical accident, though this seems unlikely, as the comparison of two different schemes in 
Tables ^ and VI show a definite stability of th ese r esults. 
For completeness, we also studied in Table VHI similar PA results, but obtained when truncating to 1/N 2 the 



perturbative coefficients beyond second (perturbative) order, and thus consistently compared to the 1/N expan- 
sion of the exact mass gap. (This is motivated from the fact that such coefficients were generated from the 1/N 2 
information [|38| exact for the RG functions, but the higher order dependence 0(1/N 3 ), artificially generated by our 
algebraic procedure (detailed in Appendix ^), is clearly not the exact one). The numerical results are very good, 
especially in the MS scheme, where it almost indicates a (numerical) convergence as the perturbative order increases. 
Though the previous numerical behaviour may be considered quite satisfactory, the choice of PA is not unique, and 
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N (0(2N) model) 




£*2 


*, 


d 4 
(2bo) 4 


2 


1.125 [0.375] 


-0.477 [0.085] 


2.16 [1.514] 


6.39 [-2.98] 


3 


0.469 [0.156] 


-0.216 [0.077] 


0.810 [0.439] 


0.508 [-0.698] 


5 


0.211 [0.07] 


-0.094 [0.045] 


0.321 [0.158] 


-0.012 [-0.211] 


8 


0.115 [0.038] 


-0.049 [0.027] 


0.162 [0.078] 


-0.04 [-0.093] 


oo 















TABLE VII: M p /A perturbative coefficients of 1/F, Eq. (^), in the scheme g 2 -> g 2 (l + bi/b g 2 ) \ and in the original 
scheme respectively [in brackets] up to fourth order. The third and fourth orders include information from the exact l/N 2 RG 
dependence. 



TV (0(2N) model) 


l/N 2 "exact" 


order 2 


order 3 


order 4 


2 


1.7293 




1.781 [1.761] 


1.778 [1.757] 


3 


1.4246 




1.438 [1.4262] 


1.434 [1.4236] 


5 


1.2252 




1.254 [1.2264] 


1.251 [1.2244] 


8 


1.1304 




1.160 [1.1322] 


1.1575 [1.1307] 


oo 


1 




~ 1 


~ 1 



TABLE VIII: M p /A at l/N 2 from a Pade approximant P[2, 3], in the scheme choice (4.8), optimized with respect to the scale 
parameter a. Only the l/N 2 expansion of the n > 3 perturbative coefficients is taken into account. The (unoptimized) MS 
values (i.e. a = 1) are indicated in brackets. 



different PA may indeed give quite different values of the mass gap 12 . For instance, when the order increases, one can 
obtain in some cases more than one extrema with respect to the scale dependence (7.E). Moreover, as already men- 



tioned the scheme is not completely fixed by our procedure, and when considering the full scale and scheme resulting 
arbitrariness, as well as all possible PA at a given perturbative order, one eventually finds numerous optima, so that 
it is difficult to decide which one is closest to the exact result without prior knowledge of the latter. This reflects 
the fact that our PA, constructed from the standard perturbative expansion, are still not able to get rid completely 
of the usual large freedom due to the scale and scheme arbitrariness in renormalizable theories. On the other hand, 
calculating directly (i.e. without optimization) in t he M S scheme (in which incidentally the exact mass gap results 24 
were obtained), we observe from Tables |v|, [VJ and VIII that the PA P 2 ,3 results, thus depending on five parameters 
appear optimal. Clearly, PA of lowest orders are inappropriate to correctly match the complete information from 
perturbative expansion beyond second order. Second, and quite interesti ngly, it happens that the higher order PA 
almost systematically have poles at u > 0, rendering the integration in Eq. (7.1) ill-defined (and therefore numerically 
unstable). In fact, as explained in Appendix B, the expected factorial behaviour at large perturbative orders is rooted 
in our PA construction, which involves at perturbative order n a contour integral of the form: 



1 

2tt7 



dv e v ln p v 



(7.6) 



P 



which is, not surprisingly, very similar to the standard renormalon behaviour Eq. ( |3.2| ) discussed in section 3. In 
other words, th oug h the PA gives a resummation method clearly different from the Borel method, the resulting 
integral in Eq. ( |7.l| ) remarkably shares with the latter some similar properties, also exhibiting the singularities of 
large perturbation orders. But, as already mentioned, the choice of PA is not uniquely fixed e.g. by purely physical 
considerations. All this gives intrinsic limitations to such PA approach, which make it inadequate to check rigorously 
numerical convergence at very high orders, as it does not take advantage of the obtained Borel convergence properties, 
discussed in previous sections. It should likely be possible to define different PA, which may better exploit the Borel 
convergence properties, avoiding in this way the usual factorial divergences, though we refrained to try such analysis 
here. It is in fact simpler to cons id er ot her kinds of numerical approximations, more directly related to the good 
convergence properties of Eq. (5.4)-( |6~16 ), as we examine next. 



12 A detailed systematic analysis of various orders of PA is performed in |39 
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C. Borel Resummation results 



Finally in this sub-section we investigate to some extent the Borel convergent resummation constructed in sections 
5 and 6. We first show in Table IX the M p / A values obtained in the rh — > limit of the (second RG order) direct 
Borel sum expr ession Eqs. (4.11), ( 4.12| ). We give here the absolute value |M P |/A due to the fact that, at second 
RG order, Eq. ( 4.1 2| ) picks up an imaginary part, ~ (— l)~ c . (This does not indicate an ambiguity of the Borel 
resummation: it is simply due to the branch cut Re[F] < ~C already present within the pure RG dependence, 
see e.g. Eq. (UU)). As one can see, the results are not very good, e xcept for the general trend that M p /A — > 1 as 
N — > oo. This is not much surprising since, as already mentioned, Eq. ( 4.12 ) is based only on the large order behaviour 
of the perturbative coefficients. We may eventually expect better results if we could include in a consistent manner 
the exact N dependence of the first few perturbative order coefficients rather than their asymptotic approximations. 
This appears qualitatively similar to the simplest oscillator case, where the direct Borel sum is known to be a poor 
approximation to the exact energy levels. 

Seco nd, w e consider the non-linear DE-VIP resummation method, explained in detail in section 6.2, with the resulting 
Eq. fl6.13|) for the large order behaviour, but now truncating this asymptotic behaviour and taking the first two exactly 
known perturbative coefficients in Eq. (2.12). Then, optimizing with respect to the new mass parameter mf, we obtain 
the results shown in table |x|. As one can see, those results, which now largely exploit the good convergence properties 
of Eqs. ( 6.15| ),( |6.16D (and the true N dependence of the first perturbative coefficients) are in more reasonable agreement 



with the exact mass gap, than the naive optimization results of section 7.1, or the Borel sum results in Table [X. But 
they are not as good as the PA results in section 7.2. As already mentioned, this is likely due to the fact that, even 
if our alternative series are formally Borel convergent, the use of such Borel summations is numerically limited. Still 
it is quite satisfactory that within this non-linear and (formally) convergent alternative expansion, there always exist 
optima, reasonably close to the exact results, which are moreover relatively flat, and stable against inclusion of higher 
perturbative orders. 



N (0(2N)) 


Mf, xact /A 


\M P \(m = 0)/A 


2 


1.8604 


3.66 


3 


1.48185 


1.69 


4 


1.3186 


1.40 


5 


1.23668 


1.27 


8 


1.1330 


1.138 


00 


1 


~ 1. 



TABLE IX: M /A values as obtained in the massless limit of the second RG order direct Borel sum expression Eq. (4.12) 



N (0(2N)) 


Mp xact /A 


M7'/A 


2 


1.8604 


1.820 


3 


1.48185 


1.375 


4 


1.3186 


1.249 


5 


1.23668 


1.187 


8 


1.1330 


1.108 


oo 


1 


~ 1. 



TABLE X: Optimized values of M p /A with respect to m' within the non-linear DE-VIP method Eq. (5.13), replacing the first 
two perturbative coefficients with their exact expressions. 
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VIII. CONCLUSION AND PROSPECTS 



In this paper, we have re-analysed in some detail a previous construction where an alternative perturbation 
expansion, based on a physically motivated self-consistent RG mass solution, can be Borel convergent in a range 
of the expansion parameter relevant for the massless limit of physical quantities in asymptotically free theories. 
The perturbative infrared renormalon ambiguities of an AFT, usually preventing unambiguous resummation of the 
standard perturbative expansion, are expected to disappear (or more precisely to cancel out with non-perturbative 
contributions) in truly non-perturbative calculations. However, such explicit cancellations are generally not possible 
to work out explicitly, except in some particular 2-D models and/or approximations, where exact non-perturbative 
results are known, e.g., at the next-to-leading 1/N orderQ, |34|, [35], [36] . In contrast, one of our main result is that 
the alternative expansion in 1/F near the relevant massless limit of the AFT can exist in two modes, thanks to the 
infrared properties of F: 

i) in the standard mode, corresponding to real F > 0, and matching the usual perturbative expansion for rh 3> A, the 
perturbative expansion alone has to be necessarily completed as usual with "non-perturbative" power corrections. 
This is illustrated explicitly with the exact 1/N calculation of the GN mass gap[p4|, as discussed in section 3. The 
net result after cancellations, Eq. ( |3.6| ), is an expression of the mass gap which depends neither on "perturbative" 
nor "non-perturbative" contributions. 

ii) In the "alternative" expansion mode, corresponding to F < in the simplest case (first RG order) where 
Aq = 1, it gives the sign-alternation in perturbative coefficients required for Borel convergence. More generally, in 
an arbitrary RS, F has complex branches near the massless limit rh — 0, so that the usual infrared renormalon 
singularities in the Borel plane are moved away from the real positive axis range of Borel integration. In this 
alternative mode the perturbative series is directly Borel summable and non- ambiguous. There is, therefore, no 
explicit non-perturbative (power correction) contributions needed in principle in this mode. We emphasize that these 
are generic results for a renormalizable AFT, encoded in the basic RG properties of the implicit function F. It should 
be stressed, however, that our construction does not solve in general the non Borel summability of AFT for arbitrary 
perturbative expansions in g(p 2 ). Rather, it is simply that the alternative expansion relevant for the mass gap, 
that is for fixed p 2 (or similar on-shell Green functions relevant to the massless limit), is performed in a neighbor- 
hood of rh — , i.e. F — or F = —C\ which is away from the usual singularities in the (complex) coupling g(p 2 )|Q, p"9| . 

In a second stage, we performed a 5-expansion (variationally improved perturbation), combined with the previous 
alternative expansion. The latter reorganization of perturbative expansions in AFT makes those particularly conve- 
nient for a 5-expansion approa ch, s ince they are much more similar to the oscillator energy levels expansion, exhibiting 



a dependence on rh v /A (Eq. ( 2.18 )) which is power- like (rather than log- like) for a sufficiently small mass rh v . This 
also explains intuitively the peculiar damping mechanism of factorial di verg ences, obtained when the trial parameter 
rh v /A is order-dependently rescaled, in analogy with similar results [|Tl[ [l2|| for quantum mechanics. Yet, unlike the 
oscillator the sole rescaling of the mass is insufficient (in the linear S expansion) to obtain a readily convergent series, 
due to the reminiscence of the RG dependence in rh v /A, making the resulting power series in rh v /A much more 
involved than the corresponding oscillator ones. (Intuitively also from pure dimensional analys is, a difference is that 
the oscillator energy level expansions in g/m 3 makes it possible to rescale the mass, cf. Eq. ( |5.3| ), such that it can 
overcompensate the perturbative factorial behaviour forpT|, |l2] ] 1/3 < 7 < 1/2, thus giving a reorganized series with 
an infinite convergence radius. While for a renormalizable D > 2 field theory, it appears only possible, at least within 
our approach, to just compensate exactly the perturbative factorial behaviour, while still being compatible with the 
massless limit.) But, when combined with Borel resummation, the usual infrared renormalon singularities at finite 
to > in the Borel plane are rejected towards to — > +00, as a consequence of the damping of renormalon factorials 
from the ^-expansion. Still, for real F > the Borel integral does not converge at t — > 00, while convergence can be 
obtained again in the range F < or complex. In this case, the ^-expansion combined with the infrared properties 
of the 1/F expansion can lead to a n im proved (fastly Borel convergent) expansion. We stress also that the linear 



DE-VIP taking the form (5.4), and (6.1) when combined with the Borel method, is only one among various similar 



resummation means. In fact the improved convergence properties do not depend on the detailed properties of the 



contour integrals here considered, e.g. Eq. (5.4) [though those have the advantage of giving rather simple and tractable 
expressions in the massless limit and for Bo rel tr ansforms Eqs. QS.lj )]: a non-linear version of the DE-VIP expansion is 



also possible, as we investigated here, Eq. ( 6.16| ), which can even lea d to a directly convergent series. More generally, 



per form ing a "brute force" ^-expansion on e.g. the mass gap Eq. ([2.6]), and rescaling the trial mass m v according 



to (5.3), replaces Eq. (5.5) and subsequent results with more complicated series, but having similar asymptotic and 
(Borel) convergence properties (see Appendix E). 

We have then performed a rather detailed numerical application of this construction for the 0(N) GN model, 
taking the mass gap (exactly know for arbitrary N from other non-perturbative approaches specific of integrable 
models) as a test of our method. Though the construction is quite general, it is restricted in practice by the lack of 
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knowledge beyond the few first orders of the genuine perturbative coefficients in most D > 2 field theories. In order 
to take into account as much as possible these uncertainties, we considered various different approximations: 
-direct numerical "PMS" optimization with respect to the trial mass of the alternative l/F expansion form, or the 
(5-expansion form of the mass gap, with the exact first two perturbative coefficients; 
-Pade approximant resummations, including higher order RG information at the 1/N 2 level; 
-"Borel- inspired" resummation, assuming ony the leading asymptotic behaviour of the perturbative coefficients. 
We obtained numerical results which are more or less consistent with each others, and showing in some cases a 
very good agreement with the exact results, at the few percent level or less, even for relatively small N values 13 . 
The best results are obtained from a certain class of Pade approximants. These PA exhibit by construction some 
of the analytical properties of the Borel sums, at least qualitatively, but can more conveniently incorporate the 
exact (non-asymptotic) lowest orders perturbative dependence as well as systematic higher order corrections. The 
optimization results are also in general better than those obtained from direct Borel sum expressions. (This is to some 
extent analogous with the Borel summable oscillator case, where the direct Borel sum is of not much practical use 
for a numerical determination of the energy levels, more accurately obtained either by Pade-Borel resummations ||, 
or the ODM [[| and related optimized delta-expansion JTl], |l2) methods). We cannot conclude from our numerical 
analysis, however, that a definite and rapid numerical convergence is obtained in field theories, in contrast with 
similar analysis performed for the oscillator energy levels. We think this may be mostly due to the large differences 
in renormalizable field theories between the true perturbative low order coefficients and their asymptotic behaviour, 
together with the fact that there are some practical limitations preventing to perform our analysis to arbitrary high 
orders, unlike the oscillator case. Nevertheless, we are confident that it should exist a better numerical exploitation 
of our construction, able to approach arbitrarily close to the exact mass gap from purely perturbative information, 
by taking advantage in a more efficient way of the formal Borel convergence properties here discussed. 
Finally, since the construction is valid a priori for any AFT, and can rely only on the perturbative information 
which is available in many models, we argue that such a summation recipe can provide a well-defined basis to 
estimate more precisely some of the XSB order parameters in more complicated, 4-D theories like QCD or other 4-D 
models with a nontrivial XSB structure. In QCD, rather than the quark pole mass which is ill-defined due to the 
confinement, the quantities of much interest are certain (gauge-invariant) condensate operators which are the order 
parameters of the (dynamically broken) chiral symmetry. These can be inferred in the chiral limit m — ► similarly 
from this alternative expansion |f6|. A detailed investigation in QCD is however beyond the scope of the present paper. 
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APPENDIX A: FORMAL RG RESUMMATION AT ARBITRARY ORDERS 



In a scheme where &j,7, 7^ for i > 2, one can generalize formally expression (2.6). We obtain 



M p {m) 



F 2 A [C + F 2 }^^ (2b F 2 )^ 



(Al) 



where by cons truction we kept for convenience a form similar to the second order, by integrating the basic RG 
equation ( [2. 13 ) separating explicitly the 60, b\, 70, 71 dependence from the higher order coefficients, that are integrated 
perturbatively, as they give a simple polynomial dependence. Thus in (Al) we explicitly separate a resummed RG 
dependence, the term F 2 A (C + F 2 )~ B , from the remnant which is systematically expanded in powers of F 2 ~ n , where 
F 2 is given by the exact second or der r elation in (2.£): the form of (2.6) as a function of F is not unique at arbitrary 
orders: the important point with (|2.6|) is that all higher order RG dependence on b n , 7n RG coefficients can be put 



Recently another variant of the PMS optimization approach, so-called "source inversion" method|4l|, also gave reasonably good numerical 
results for the GN model. This latter approach, which shares some similarities with our method, does not address, however, the important 
problem of the large perturbative order behaviour. 
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into the form of contributions to the coefficients d n of F 2 n , as explicitly shown below. Accordingly we can write 

d n = c n + r n (A2) 

to distinguish contributions from the true perturbative part, c„ (i.e. the non-log part of the n-loop perturbative 
graph) from r n , the contribution for n > 2 originating from higher order RG-dependence (in a scheme different from 
Eq. (2. 2), (2. 3)) due to the re-expansion in powers of 1/F implicit in (2.6) or (Al). 



At arbitrary order n, F(rh), A and the invariant mass rh are different from the corresponding second RG order 
expressions. Their expressions can be again formally derived order by order in perturbation. For instance 



F(m/ A) = F 2 + <*n/F2 + £ 5 n /F% 



(A3) 



n>l 



n>2 



m = m 2 cxp[-— 2^ o n g \ 

n>2 



(A4) 



where 



b 2 hb 2 



26g Abl ' 



5 2 = 



(A5) 



and m 2 in Eq. (A4) designates the second RG order scale invariant mass, defined in Eq. ( |2.7| ). Similar ly, A at arbitrary 
RG orders n involves perturbative corrections with respect to the (universal scheme) expression in (2.10): 



A = (j. exp(- 



)ib o9 \nT c [i + JV] c [i + EW>)] 



2b g 2 (fi)' 



(A6) 



with 



Ai 



A, 



2b 2 b 3 - 4boMg - b\ 



(A7) 



All the above relations expresses consistently the fact that A is scale invariant and rh is scale and scheme invariant at 
a given n order. Moreover, the connection between the exact second order part F 2 of F, which is not perturbatively 
expanded, and the exact second order part rh 2 of rh, is such that one still has the property 



m,F- A {C + F)- 



const. 2-° A + )] 



(A8) 



when rh — > 0, where A = ji/(2b\), thus maintaining the infrared properties of the massless limit discussed at first 
and second RG orders in section 2. 



APPENDIX B: OSCILLATOR ENERGY LEVELS AND BOREL SUMM ABILITY 

In this appendix we briefly elaborate on a perhaps unusual picture of the anharmonic oscillator energy levels, in 
order to illustrate in a simpler case the existence of two different perturbative expansion modes, with one mode 
needing additional non-pcrturbative corrections while the other mode is directly Borel summable. 
We recall that the oscillator is described by a D = 1, g<f> 4 massive scalar field theory |fi"3"|: 

A«, = \{d t 4>? - - ^ (Bi) 

The energy levels have a perturbative expansion 1 



We shall only consider the ground state energy level to simplify, which is sufficient for our purpose. 
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where the coefficients can be calculated to arbitrary orders, and have the well-known asymptotic behaviourjll 



_ ( _i)" ( ^)V2 ( ! rr[n+1/2] (i + ^I)) 

7r J 4 n 



(B3) 



Because of the sign- alternation of the coefficients in Eq. (B3), the series is Borel-summable[|5|. Consequently, there 
are no additional "non-perturbative" contributions, and this series can be uniquely Borel-resummcd to represent the 
(real part of) the oscillator ground-state energy. Explicitly, the Borel sum reads: 



E (g) 



(6) 



1/2 



.9 



dte 



-1/2 



(B4) 



Notice that the result Eq. (B3), originally obtained from a WKB calculation |13|, was in fact later shown to be derivable 
from an instanton-based calculation, the vacuum being unstable for g < with a corresponding tunneling process. 
More precisely, a classical calculation gives |i7j for the instanton contribution to the imaginary part of the ground-state 
energy: 
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4 4| m | 3 

e 3 s 
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1/2 



-<J 



(B5) 



where g < 0. 

Now let us assume that our only knowledge would consist of the purely perturbative information, namely the expansion 
in Eq. (B2), and consider formally changing the sign of the coupling there: g — > —g. Obviously this cancels the sign- 
alternation of the coefficients, so that the corresponding new series is no longer Borel summable: more precisely 
instead of Eq. (B4) one obtains an integral 



E (g) 
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1/2 



1.91 
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V 1/2 

4 ; 



(B6) 



which is ill-defined due to the cut at t = +4/3 on the integration range. Nevertheless, we can evaluate the ambiguity 
that this implies: by defining the ambiguity by the (half) difference of the two possible contours avoiding the cut by 
above (resp. by below), a straightforward contour calculation gives an ambiguity: 
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4i 

7^ 
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1/2 



(B7) 



which accordingly is to be interpreted as an additional non-perturbative contribution needed in this unco nven tional 
picture. Not suprisingly, one recovers in fact consistently the non-pertur bati ve instanton contribution Eq. (B5), with 
g < 0. (NB we r ecal l that actually, it was the instanton solution Eq. ( |f35| ) that lead to derwe|37|] the asymptotic 
behaviour in Eq. (B3), so that Eq. (B7) is nothing but a consistency check, though relevant to our argument.) Thus, 
the two perturbative expansion modes, the directly Borel summable one with g > or g < respectively, can be both 
consistent, provided one correctly identifies the necessary additional non-perturbative contributions in the g < case. 



APPENDIX C: PADE APPROXIMANTS AND CHIRAL LIMIT 



In this appendix we give some technical details on the construction of the Pade approximants (PA) used in sec. VII B . 
Similar PA were used in ref. H]. As mentioned, the aim is first to have a we ll-defined (finite) massless limit, F — > 0, 
to/ A — > , wh en including the purely perturbative expansion series. Eq. ( |2.8[ ) suggests that for the contour integrand 
e.g. Eq. (^4), where F = F[m" v], one reintroduces instead of m" = fh/A an (inverse) "effective coupling" variable 
u defined by 



u A e u (C + u) 



B-C 



(CI) 



In the usual perturbative regime, m" 
asymptotic expansion for u — > oo: 



oo, in terms of it, the In In. . . lnm" disentangle, and Eq. (|2.8|) admits an 



F[m" v] ~ u 



\nv_ _ (A + B-C) ^ F _ {B-C) (1 + C/F) . 

u u u u (1 + c/u) ' 
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(A + B-C) lm 



+ 0(u~ 3 ) 



(C2) 
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p 


I p 


1% (contour) 





1 





1 


-7s ~ -0.577216 


1 


2 


7 | + C[2] ~ 1.97811 


1.15443 


3 


-7s - f 7b - 2C[3] ~ -5.44487 


3.93527 


4 


7b + ^ 2 Tl + + 8 7bC[3] ^ 23.5615 


-1.00806 


5 


■ ■ ■ ~ -117.839 


-19.9846 



TABLE XI: Results of integrals Eqs. (C5), (C4) for the first few p values 



where the higher order terms are polynomial in lnu, and evidently depends on the RG coefficients A,B,C. Using 
Eqs. (CI), (C2) within the contour integral Eq. (5.4), around the cut at Re[v] < defines the mass gap as an 
asymptotic expansion: 
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A 



const. e u u~ (c+1) x 
1 



A Q + -{{A + d t + d x A - A (A -C) + ~f E A {A + 1)) + 0{ u - 2 ) 
u 



(C3) 



where Aq = A + B 1 and di are the original perturbative coefficients of Eq. (2.6) in a given scheme, including eventually 
the 1/N 2 RG-dependence (see Appendix O). The overall constant includes any constants from e.g. RG dependence 
(such as typically the 2~ c factor in Eq. ( |5.4| )), or depending eventually on the choice of renormalization scheme. In 
order to derive perturbative expansion such as Eq. (C3) at arbitrary higher orders u~ p , M p /A has to be systematically 
expanded using 



2iri I 7r ' J 

J 3=1 



p-2 



(C4) 



where (?) are the binomial coefficients and 



In = 



(C5) 



can be evaluated exactly in terms of Euler's constant ~/e an d the Riemann zcta functions £[n]. We give in table XI the 
first few p orders for those integrals. Note that the contour integral results 7£ in (C4) are real. We remark that the 
basic integral I p behaves as ~ (— l) p pi as p increases, and similarly the final contour integral F% behaves as ~ (p— 1)!, 
with no definite sign. 

Taking now the derivative with respect to u of the logarithm of this expression for M p /A yields the final power series 
in 1/u 



du 



C + l 



0{u 



(C6) 



suitable for PA analysis. Similarly to jL4|, we thus consider PA P[ p , q ] as rational fractions of polynomials: 



(C7) 
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with coefficients s n , t n , uniquely fixed fro m requiring that the perturbative expansion, for u — > oo, matches the 
perturbation series obtained from Eq. (|C6|), and that the massless limit u — > is finite. Actually since our PA are 
constructed from d ln[M p /A]/<9u, one further step is needed to recover the mass gap itself. We get more precisely 



M p M p 
In — - — (u — > 0) = In — - — (u — > oo) 
A A 



du P [p , q] (u) 



(C8) 



where the initial condition 



ln-^-(u -» oo) = ln[const.A ] + u - (C + 1) lnu + 0(l/u) 



(C9) 



is uniquely fixed because it is entirely determined by the first order of perturbation theory. Its (perturbative) diver- 
gences for u — > oo are exactly compensated by the corresponding ones from the first terms of the Pade approximants, 
and we finally obtain 



M p 

— = const. Aq exp 



1 - 



du P [Pi9 ] (u) du [P\p,g] (u) - 1 + 



C+l. 



(CIO) 



which is used in the numerical analysis in section VII B 



APPENDIX D: FINITE 1/N 2 CORRECTIONS FROM SCHEME CHANGE 

In this appendix we give relevant formulas for the RS change which defines the (approximated) 1/N 2 corrections to 
the finite perturb ative coefficients of arbitrary orders. These are used for the numerical analysis of Pade approximants 
in section VII B 

The first step is to extract from the results of ref. |38| the perturbative expansions of the RG functions /3(g) and 7 TO (<?) 
to arbitrary orders in the MS scheme, which serves as the basis of our RS change. In fjsl, by studying the theory 
near the (perturbative) fixed point g c of the e-dependent beta function: 



P(g 1 e)^e + J2(-2b l )(g 2 ) 



2+i 
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(Dl) 



were derived the exact 1/N 2 expressions, as function of e = D — 2: 

/3'(g c )^^(g c ) = -2^-1 -(2^-1)^ 
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N 2 



(D2) 



and 



lm(g c 



D-2 



m + xi . V2 + X2 



N 



N 2 



(D3) 



In (D2) and (D3) /i = (D + e)/2 and the coefficients rji, Xi and Xi are given analytically in terms of the Gamma 
function T[z] and its first two derivatives, ^[z] and \& [z\. We refer to Ref. for details on this derivation and 
explicit expressions of rji, Xi, Xi which we essentially follow except for slight differences of conventions e.g. in the 
definition (Dl). 

The first three-loop orders hi, ji up to i = 2 are known for arbitrary N in the MS scheme |2S)|. Next, by simply 
comparing the expansion of (3 (g c ) to arbitrary orders in e, using Eq. (Dl): 
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(where b n appears first at order e n+1 ) with the same e expansion obtained from the exact 1/A^ 2 expressions (p2|), 
(p^), which gives e.g. for (3'(g c )' 
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it is straightforward to obtain e.g. for the n > 4 loops unknown RG coefficients: 
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and 
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where we refrain to give higher orders which involve quite lengthy expressions 15 . Ne xt step is to incorporate this exact 
l/N 2 information as corrections to the 't Hooft scheme mass gap expression ( |2.6| ). To this ai m, w e exploit the fact 
that the specific change of scheme implie d by pass ing from the MS RG functions Eqs.(Dl)-(D7) to the truncated 
ones in the two-loop 't Hooft scheme Eqs. (2.2),( |2~3| ) can be achieved by appropriate redefinitions of the coupling g(fj.) 
(equivalently redefining 1/F) and redefinitions of the mass m(/i). The most general perturbative such RS change (at 
fixed scale [i): 
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m — > fa = m Z m (g) — m(l + Big 2 + B 2 g A + •••) 



■0(9) ^ 



implies fromffO| 

7m (5) = 7m (g) ~ 25 /%) afs \nZ m (g) 
relations between the RG coefficients in the two schemes [Ell 



(D8) 



(D9) 



7i 
k 
72 



71 + 2b B 1 - 

b 2 - Aybi +b (A 2 -A 2 ) 

72 + 2biBi + 2bo(2B 2 - B\) 



7 A 2 



dx =d 1 -B 1 , d 2 = d 2 - di(Ai + Si) - (B 2 - 70.B1) 
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and corresponding changes in the purely perturbative coefficients of (E76J) : 



(D10) 



(Dfl) 



The RS transformation is uniquely fixed, apart from 71, by requiring the new scheme to be the above l/N 2 MS 
Eqs (D6), (D7) and the old scheme to be the two- loop truncated RG coefficients. There is however one missing 
piece in this construction: the exact 1/iV 2 dependence above is the one of the RG coefficients, while the genuine 
(non-RG) finite parts, are actually unknown for the GN model beyond two-loops. Thus, in addition to fix uniquely 
the perturbative contributions d n in Eq^^) in the new scheme, we assume that these are negligible beyond two 

for n > 2. Whether this assumption is a good approximation or not 



in Eq.(|2 

loops, in the original ('t Hooft) scheme, d n 



can only be decided by the numerical analysis, see section 7B. 



Actually, the brute force RS change as implied by Eqs.(DlC), (Dfl) becomes rapidly algebraically involved as the 
expansion order increases. A more efficient method can be devised]39||, more convenient to study Pade approximants 
of relatively high orders. Noting that the RS transformation is uniquely fixed, up to 71, by A n and B n in Eqs.([D§|), 
one can interpret the RS transformation as a one-parameter r variation in the RS parameter space, with "boundary" 
conditions A n {r) and B n (r) for the initial and final scheme, r = and r = f, respectively. Using that the bare 
coupling and mass are RS invariants, one defines 



dg_ 

dT 

dm 



(D12) 



All algebraic calculations were performed with Mathematicap2[ . Note in Eqs. (|i5|)-(|D7|) that N now refers to the O(N) model 
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with perturbative expansion 



a(g,r) = -a g 3 - a x g 5 
S m {g,r) = D g 2 + D ig 4 - 



(D13) 



similarly to the above ordinary RG functions (3(g) and J m (g)- While the latter parameterize the variation with the 
scale /j,, the functions in Eqs. ( D12| ) parameterize the variation in the most general RS parameter space (which is 
2n + 1 dimensional at perturbative order n) , with r playing a role similar to In /x, and with 2n boundary conditions 
fixed by the initial and final specifi ed sc he mes. This is completely equivalent to the brute force above RS change, 
but more convenient, because Eqs ( D12| ), ( D13 ) can be resummed to some extent, for the leading and subleading 
dependence, similarly to the well-known RG resummations properties of (3(g) and 7 m (<?)- More precisely, imposing 
for simplicity that the coefficients ai and Di do not depend on r, A n {r) and B n (r) are series expansion in t and we 
also have the following useful relations between the RS change functions and the ordinary RG functions: 



dT 



da 
d\i 



d'Jra 
dT 



d5 m 



(D14) 



which gives by perturbative expansion 



-r-(b r -i) = 2^2(2n - r)a r _„_i6„_i(r) 



(D15) 



n=l 



and 



dT 



(7o 





r-l 



dT 



(7, 



r-l ) 



^2 2 ^[ar-n-l7n-l(T) - 



(D16) 



n=l 



Integration of these equations leads to appropriate recurrence relations ]39| : 

r-l „ T 

6 r _i(r) =2V(2ti-r) a r _ n _i / ds6„_i(s) 

n=l Ja 

7r-i(r) =Y^2ri a r _„_i / ds7„_i(s) - D„_i / ds7„_i6 r _„_i(s) 
defining order by order all the needed RS information. 



(D17) 
(D18) 



APPENDIX E: STANDARD DELTA-EXPANSION 



For completeness we briefly examine here the DE-VIP convergence properties in our framework by applying to 
the mass gap expression Eq. (2.6) the more standardfl9|-]l2]] order by order ^-expansion without our specific contour 
integral resummation prescription illustrated in sections 5 and 6. To this aim we first define the substitution 



rn 

,2^ , X „2 



m v (1 + 8 -) = mX(l - S /3{X)) 



where X = m v /m and (3{X) = 1 — l/X, to apply to 



(El) 
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expanded to order S N . (Without much loose of generality we again took for simplicity in Eq. (E2) the first RG order 
expression of the resummed RG dependence). After straightforward algebra it gives 



M p /A 



47T b 



N N-qN-q-p 

EE E 

q—l p—0 r—0 



T\p + q]{p + q + A){q + Af- 1 T[r + q/A] 



APT[l+p] 



T[r + l}T[q/A 



) (m"X)-i/ A (P(X)Y 



(E3) 



The basic calculation is similar to the one that lead to Eq. (5.5), except for the factor l/T[l + q/A] less (which in ( p.q ) 
originates from the contour integral) , and the third summation on r in addition, which corresponds to the expansion 
of [X(l-S /3pO)]-«/ A . Note also that the delta -expansion is truncated at order TV, and 5 — ► 1: since the perturbative 
term g 2n — (g 2 ) p+q — + S p+q (g 2 ) p+q , it explains the third summation upper bound N — q — p. 
Note furthermore that here the expansion of F around is rigorously valid, as long as 



mX(l - S P(X)) < e - A A A 
(where the right-hand side corresponds to the convergence radius of the power expansion form of F). 



(E4) 



It is clear that the terms of the third series, though considerably complicating t he a lgebra, do not play any role in 
the asymptotic/convergence properties. Looking at the asymptotic behaviour of (E3) it is easily seen, by rescaling 

X = m v /m as 



X -> iV 7 X . 



(E5) 



that the series behaves like the one in (5.5), for 16 A < 7. 

The series can next be "Borelized" in a way similar to what is described in section 5 above. It gives more complicated 
Borel integrals, but having similar asymptotic and Borel convergence properties, in particul ar for Re[m v ] < 0. We 
refrained however to attempt any numerical analysis based on the complicated series in Eq. (E3), which appears less 
convenient for this purpose than the construction mainly discussed in the rest of the paper. 
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